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Preface 


The treatment in this book emphasizes at all times that a vector is an 
entity in itself and should not be considered merely as a triple of 
numbers. Consequently, proofs by means of Cartesian coordinates 
are avoided. 

The first chapter contains the basic elements of vector algebra, 
which are applied in the following chapter to the geometry of the 
straight line and the plane. The third chapter discusses the funda¬ 
mentals of vector calculus, and applications to differential geometry 
are made in the fourth chapter. The contents of this chapter, with the 
exception of the Frenet formulae, curvature and torsion, are required 
for later developments of the theory. 

Line integrals, surface integrals and volume integrals play an 
important role in vector analysis. Many students will meet these 
concepts before they have had formal courses in them. Consequently 
the fifth chapter presents a concise account of these topics. 

The following three chapters are respectively devoted to the 
gradient, divergence and curl of a vector. Of these, the first is defined 
in terms of the directional derivative whilst the latter two are defined 
by limits of integrals. In this way the theory is developed inde¬ 
pendently of any coordinate system. The operator V is not introduced, 
as the author has observed that this is a dangerous tool in the hands of 
the inexperienced. 

There follow chapters on Stokes’s theorem, Green’s theorems and 
orthogonal curvilinear coordinate systems. 

The concluding chapter provides a link between vector and tensor 
analysis. It includes a discussion on the physical components of a 
vector as distinct from the components referred to a general basis. 
There is some confusion in the literature of vectors about the validity 
of the equation curl curl a = grad div a—da referred to a curvilinear 
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coordinate system. Some authors in fact define Aa by means of this 
equation. However, at this point a so-called natural basis is introduced 
with respect to which this equation can be deduced. 

No attempt has been made to apply vectors to problems in mathe¬ 
matical physics. Indeed, a comprehensive treatment which would do 
justice to those applications would necessitate a threefold expansion 
of this book. 

A representative selection of examples is inserted at appropriate 
places and answers are provided at the end of the book. 

During the preparation of this book, I have been greatly indebted 
to Dr. J. H. Wilkinson who read the manuscript and made many 
valuable suggestions which have been incorporated in the work. 

My thanks are also due to my colleagues Dr. M. G. Smith and Mrs. 
R. J. Church, the former for a critical reading of the manuscript and 
the latter for the supply of a number of examples. For help in reading 
the proofs I wish to thank my wife and elder son. 

Finally, I wish to express my appreciation to the staff of the 
publishers for their cooperation and for their help with the diagrams. 

B. S. 
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CHAPTER 1 


Vectors 


1 Vectors 

Physical quantities such as mass, temperature and work are measured 
by numbers referred to some chosen unit. These numbers are called 
scalars. 

Other quantities exist such as displacement, velocity, acceleration 
and force, which require for their complete specification a direction 
as well as a scalar. These quantities are called vectors and may be 
represented by a straight line with an arrow. 

p --- Q 


Fig. 1 

Formally, a vector is a directed line segment. The vector depicted 
in Fig. 1 from the point P, called the initial point, to the point Q 

will be denoted by PQ . It has the length or magnitude PQ and the 
direction from P to Q as indicated by the arrow. 

In the special case when Q coincides with P, we refer to the zero 

vector PP denoted by 0 . This vector has zero magnitude but indeter¬ 
minate direction. 

In addition to the notation PQ for a vector, it will be convenient 
to designate a vector by a letter in bold type such as A, a, a, etc. 

The magnitude of a vector a will be denoted by either |a| or a and 
it is well to emphasize that the magnitude of a non-zero vector is a 
positive quantity. 

A vector of unit magnitude is called a unit vector. The letters i, j, 
k and I, J, K are reserved for unit vectors. 

Two vectors are said to be equal if they have the same magnitude 
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and direction. Geometrically, all vectors obtained from the vector PQ 
by a translation in space are equal. 


Exercise 

1. Which of the following are scalars and which are vectors? (a) volume, 
( b ) energy, (c) momentum, (d) temperature, (e) work, (/) electric field 
intensity. 


2 Addition of Vectors 

The addition of two vectors is defined by the so-called triangle or 
parallelogram law. 



Let PQRS in Fig. 2 be a parallelogram. Formally we define the sum 

of the vectors PQ and QR to be the vector PR represented by the 
third side of the triangle PQR. 

Since the vectors QR and PS are equal, we may also define the 


sum of the vectors PQ and PS, which have the same initial point P, 

to be the vector PR represented by the diagonal through P of the 
parallelogram formed by PQ and PS. 

Let us write 


PQ = SR = a, 


PS = QR = 


b. 


By the triangle law of addition, both the sum of PQ and QR and the 
sum of PS and SR are equal to PR. It follows that 

a+b = b + a. 


That is, vectors satisfy the commutative law of addition. 
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With reference to Fig. 3 let us write 

PS = a, SR = b, RQ = c. 

By the triangle law of addition we have 

PR = a+b, SQ = b + c. 

Further applications of the triangle law yield 

(a+b) + c = PR+RQ = PQ = PS+SQ = a + (b + c). 

That is, vectors satisfy the associative law of addition and so there is 
no ambiguity if either expression is written without brackets in the 
form a + b + c. 

Figures 4 (a), (b) and (c) show that in the case when PQ and QR 
are parallel, the triangle law yields the result 

]PQ + QR = PR 

whether the arrows along PQ and QR are in the same or opposite 
directions. To avoid confusion the arrows are not drawn in the 
figures. 



Fig. 4 
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At this stage it is well to point out that a physical quantity may be 
specified by a direction and magnitude and yet not be a vector because 
it does not obey the vector law of addition. Finite rotation of a rigid 
body is such a quantity. 


3 Subtraction of Vectors 

If the vector a is represented by the directed displacement PQ, we 
define — a to be the vector represented by QP. Hence 

a —a = 0. 

Now we may define the difference a—b of two vectors a and b to 
be the sum of the vectors a and — b. That is, 

a —b = a + ( —b). 



Let PQ and PS in Fig. 5 represent the vectors a and b respectively. 
Complete the parallelogram PQRS. It follows from the triangle law 
of addition and the definitions of the negative of a vector and sub¬ 
traction that 


a-b = PQ-PS 
= PQ+lsP 
= SP+PQ 

= SQ- 


That is, the difference of two vectors can be represented by a diagonal 
of a parallelogram. 


Exercises 


1 . Show that |a+b| < |a| + [b|. 

2 . Show that [a—b| > ||a| — |b||. 
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3 . Under what conditions do we obtain equality signs in the inequalities 
of Exercises 1 and 2 ? 


4 Multiplication of a Vector by a Scalar 

If A is a scalar and a is a vector, Aa is defined to be the vector with 
the same direction as a if A is positive, the opposite direction to a if 
A is negative, and with magnitude |A| times the magnitude of a. 
Consequently we have 

|Aa| = |A| |a|. 

Note carefully that ]a| denotes the magnitude of the vector a whilst 
|A| denotes +A or — A according as A is positive or negative. Also we 
see that a/a is the unit vector in the direction of a. 

It follows immediately that 


and 


A(yua) = yx(Aa) = A fiB. 
(X+fj)a. = Xa+fxa. 


s 



Let OPQ (Fig. 6) be the triangle whose sides OP, PQ and OQ 
represent the vectors a, b and a + b respectively. Let RS , parallel to 
PQ, cut OP and OQ in R and S. Then we have 

OR RS _ OS 
OP~ PQ~ OQ 

Denoting these equal ratios by A, it follows that OR, RS and OS 
represent the vectors Aa, Ab and A(a+b) respectively. Hence by the 
triangle law of addition we have 

A(a+b) = Aa + Ab 


6 VECTORS [Ch. 

and so we are permitted to remove brackets in the usual algebraical 
manner. 

Two vectors a and b are co-directional if a = Ab, where A is 
positive; that is, if a and b are parallel and have the same direction. 

Exercises 

1. If a and b are vectors whose directions are neither parallel nor 
coincident, show that the relation Aa+/xb=0 implies that both A and [i are 
zero. 

2. If P is any point and D, E, F are the mid-points of the sides BC, CA 
and AB respectively of the triangle ABC , show that 


PA+PB+PC = PD+PE+PF. 


5 Point of Division 

Let R in Fig. 7 divide the join of A and B in the ratio AR/RB= A. 
Then AR = XRB. 



Further, let A, B and R be given by the vectors a, b and r referred 
to the fixed initial point O. Then we have 

r—a = A(b — r), 

from which it follows that 

a+Ab 

r ~ TTT 

In particular, we see that the mid-point of AB is given by the vector 
i(a+b). 

When the three points A, B and R are collinear, the previous result 
shows that a number A exists such that 

a+Ab —(1 + A)r = 0. 

Hence three non-zero numbers a, and y exist such that 

aa+jSb+yr = 0 

a+/? + y = 0. 


and 
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Conversely, given that aa+j8b + yr = 0, where a, j8 and y are three 
non-zero numbers satisfying a+j8 + y =0, there is no loss in generality 
if we assume that a +j8 # 0. Hence we have 

_ aa+j8b 

cc + j8 

and so R lies on AB and divides it in the ratio j8/a. 

The vectors a, b, c,..., 1 are said to be linearly dependent if we 

can find a set of scalars A, /n, v f _, />, not all zero, such that 

Aa+/xb+vc + • • • + pi — 0. 

Otherwise they are linearly independent. 


A 



Example 1. Prove that the medians of a triangle meet at a point which 
divides each median in the ratio 2:1. 

Let A , B and C be given by the vectors a, b and c. Then the mid¬ 
points D , E and F (Fig. 8) are given by ^(b + c), ^(c + a) and ^(a + b). 
Let G be the point which divides AD in the ratio 2:1. Then G is given 
by the vector ^(a + b + c). Similarly G lies on BE and CF and so the 
result is established. 


A 



2 —V.A. 


Fig. 9 
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Example 2. Prove that the line joining the mid-points of two sides of 
a triangle is parallel to the third side and is equal in length to half that 
of the third side. 

In Fig. 9 let A, B and C be given by the vectors a, b and c. Hence 
the mid-points E and F are given by ^(c + a) and ^(a + b). Accordingly 

FE=%(c + a) — ^(a + b) = ^(c — b). But BC= c — b and so the result 
is established. 

Exercises 

1. Show that the mid-points of the sides of a skew quadrilateral form 
the vertices of a parallelogram. 

2. Show that four points A, B, C, D, no three of which are collinear, 
given respectively by the vectors a, b, c and d are coplanar if and only if 
four non-zero numbers a, j8, y and 8 exist such that 

aa+j8b+yc + 8d = 0 and a+j3+y+8 = 0. 

3. Show that the lines joining the mid-points of the opposite edges of a 
tetrahedron are concurrent and bisect each other at the point of concurrency. 

4. The centroid G of the points A lf A 2 ,..., A n with associated numbers 
m lf m 2 , ..., m n is defined by the equation 

OG 2 m r = 2 ( m r OA r ). 

r = l r«= 1 

Show that the position of G is independent of the choice of the initial 
point O. 

6 Components of a Vector 

Let OP (Fig. 10) represent the vector a and OU the unit vector u 
which makes an angle a with a. Let the perpendicular from P to OU 
intersect it at L. The length OL, multiplied by plus or minus one 



according as OL is in the same or the opposite direction to OU> is 
called the component of a in the direction u. In either case the 
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component is a cos a. Further, the vector a cos a u is called the 
projected vector of a in the direction u. 

The projected vector of a on a plane is defined to be the projected 
vector of a in the direction of the line of intersection of the given 
plane and the plane through a perpendicular to the given plane. 

Consider (Fig. 11) the vectors OP= a and PQ = b and the unit 

vector OU=u. Let PL and QM be perpendicular to the direction 
OU. The components of the vectors a and b in the direction u are 
OL and LM respectively. The sum of a and b is represented by the 

vector OQ which has the component OM in the direction u. 


p 



Accordingly the sum of the components of two vectors in a given 
direction is equal to the component of the sum of the two vectors in 
that direction. 

This result clearly extends to any number of vectors. That is, the 
sum of the components of any number of vectors in a given direction 
is equal to the component of their sum in that direction. 


Exercises 

1. ABCDEF is a regular hexagon whose sides are of length 5 units. If 
AB =a and BC=b, find the following vectors in terms of a and b: (i) CD, 
(ii) DE, (iii) EF, (iv) FA, (v) DA, (vi) EB. Further, obtain the components 

of FA and FE along BD and AB respectively. 

2. ABCD is a regular tetrahedron whose edges are of length 3 units. If 

AB= a, DC=p and CD= y, find the vectors AD and CA in terms of a, (3 
and y. Further, calculate the projected vector of AD along CA. 
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7 Fundamental System of Vectors 

Introduce (Fig. 12) a right-handed system of mutually orthogonal 
coordinate axes OX , OY and OZ. That is, a right-handed corkscrew 
rotating through a right-angle from OY to OZ advances in the 
positive direction OX. Similarly, if the corkscrew rotates through a 
right-angle from OZ to OX or OX to OY it advances in the directions 
of OY and OZ respectively. A left-handed system is obtained if the 
direction of any one of the axes is reversed. We shall always select the 
axes to form a right-handed system. 

The three unit vectors i, j and k in the positive directions of the 
coordinate axes are called a fundamental system of vectors. 



Let a represent the vector OP. Complete the rectangular parallel¬ 
epiped with diagonal OP and edges through O along the axes. It is 
clear from the figure that OL, OM and ON are the components of 
a in the directions of i, j and k respectively and that the projected 

vectors are OL, OM and ON. 

By the repeated application of the triangle law of addition of vectors 
we have 


a = OP = OV+ VP = OL + LV+ VP 
= OL + OM+ON 


and so 


a = a*i + a, j + tf fc k, 


where a { , a, and a k denote respectively the components of a relative 
to the fundamental system. 
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Further, we have 


OP 2 = OV 2 +VP 2 = OL 2 + OM 2 + ON 2 , 

and so 

a 2 = a? + aj + a%. 


Exercises 

1. If a=4i+j—k, b=3i—2j + 2k and c= — i — 2j+k, calculate (i) 
a+b+c, (ii) |a—b—c|, (iii) a unit vector parallel to 2a—b—c but in the 
opposite direction. 

2 . Prove that the vectors i—k, — i+j + 2k and i—j —3k can form the 
sides of a triangle. 


8 Scalar Product 

The scalar product of two vectors a and b is defined to be the 
scalar ab cos 6, where 6 is the angle between the vectors a and b. 
The scalar product corresponding to the vectors a and b in Fig. 




Fig. 13 


13(a) is positive since 6 is acute, while the scalar product in Fig. 13(6) 
is negative since 6 lies between 7r/2 and it. 

We agree to denote the scalar product by a.b and so it is sometimes 
referred to as the ‘dot’ product. (Other notations may be met, such 
as a|b, (a, b).) We have 

a.b = ab cos 6. 

It follows that 

a.b = b.a. 

That is, scalar products satisfy the commutative law of multiplication. 

From the definition, a.b = 0 if either «=0, 6 = 0 or cos 6=0. In 
the latter case the vectors a and b are mutually orthogonal. Con¬ 
versely, the scalar product of two orthogonal vectors is always zero. 




12 


VECTORS 


[Ch. 

Again, we have from the definition that the scalar product a.b of 
two co-directional vectors a and b is ab. Hence a. a = a 2 . At times it is 
convenient to write a. a = a 2 . 

The fundamental system of vectors i, j and k satisfies the relations 


i i = j j =k.k= 1,1 
j.k = k.i = i.j =0./ 


( 8 . 1 ) 


If A is a scalar, we see from the definition that 


A(a.b) = (Aa).b = a.(Ab) 

and so we can write these equivalent expressions in the form Aa.b. 

Another consequence of the definition of the scalar product is that 
the component of the vector a in the direction of the unit vector u is 
a.u. 

In Section 6 we proved that the sum of the components of two 
vectors in a given direction is equal to the component of the sum of the 
two vectors in that direction. It follows (Fig. 11) that 

a.u+b.u = (a+b).u. 

On multiplication by the scalar A, we have 

a.c + b.c = (a + b).c 

where c = Au. That is, we may remove brackets in scalar products as 
in ordinary algebraic multiplication. In other words, scalar multi¬ 
plication is distributive with respect to addition. 

Let the vectors a and b have components a iy a jt a k and b if b jy b k 
respectively with respect to a fundamental system of vectors i, j and 
k. Then 

a.b = {aj. + a J + «fck) • (fyi + brf + b k k). 

Removal of brackets and use of the scalar product properties (8.1) of 
the vectors of the fundamental system yield 


a.b = afbi + ajbf+ajjbfr 


Example 1. Show that the three altitudes of a triangle are concurrent. 

Let A y B and C be given by the vectors a, b and c referred to some 
initial point. Further, let (Fig. 14) the perpendicular BY from B to 
CA intersect the perpendicular AX from A to BC at the point H 
given by the vector h. 

The conditions of orthogonality are 

(b — c).(h—a) = 0 and (c —a).(h—b) = 0. 
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From the identity 

(b—c).(h—a) + (c —a).(h—b) + (a—b).(h —c) = 0 
we have 

(a—b).(h—c) = 0. 

That is, CH is perpendicular to AB as required. 


c 



Example 2. Obtain the acute angle between two diagonals of a cube. 

Let the rectangular parallelepiped in Fig. 12 be a cube with edge of 
length unity and let the perpendicular from P to the XY plane cut it 

at V. Then the vector OP= i+j + k whilst the vector NV = i+j—k. 
Hence the angle a between the two diagonals is given by 

|i+j+k| |i+j-k| cos a = (i+j+k).(i+j-k), 

from which cos a= 1/3 and so a = 70° 32'. 

Exercises 

1. If a.b=a.c, show that either a=0, b = c or a is orthogonal to b—c. 

2. Given a=i+2j — 3k and b=3i—j+2k, (i) show that the vectors 
a+b and a—b are mutually orthogonal, (ii) calculate the acute angle 
between the vectors 2a+b and a+2b, (iii) obtain a unit vector orthogonal 
to both a and b. 

3. Show that (a.b) 2 ^a 2 b 2 and deduce that 

(«A+ ajbj + a k b k ) 2 < (af+af + a%)(bf+bf+b%). 

4. Prove that a=(a. i)i+(a. j)j + (a. k)k. 

5. Show that the perpendicular bisectors of the sides of a triangle are 
concurrent. (Hint: use the identity 

(b-c).{l-i(b+c)}+(c-a).{l-i(c+a)}+(a-b).{l-i(a+b)} = 0.) 
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9 Vector Product 

Let n be the unit vector orthogonal to both a and b and such that 
a, b and n form a right-handed system. We define the vector product 
of two vectors a and b inclined at an angle 6 to one another, where 
0 ^ 6 4: TTy to be the vector ab sin 6 n. 

We agree to denote the vector product by a x b and so it is some¬ 
times referred to as the ‘cross’ product. (Other notations in use are 

a Ab, [a, b] and a b.) We have 

a x b = ab sin 6 n. 

Since b, a and — n form a right-handed system, it follows that 
b x a = — ab sin 6 n 

and so 

axb = —bxa. 

That is, vector products satisfy an anti-commutative law of multi¬ 
plication. 

If axb = 0, then either <2 = 0, b = 0 or sin 6=0. In the latter case 
6 = 0 or 7 t and the two vectors are parallel and have the same or 
opposite senses. Conversely, the vector product of two parallel vectors 
is the zero vector. As a special case we have a x a=0. 

From the definition we have that if a.a = b.b=l and a.b = 0, the 
vectors a, b and axb form a right-handed system of mutually 
orthogonal unit vectors. 

The fundamental system of vectors i, j and k satisfies the relations 


(9.1) 


ixi=jxj = k x k = 0, 
jxk = — kxj = i, 
kxi= —ixk=j, 
i x j = -j x i = k. 

If A is a scalar, we see from the definition that 

A(axb) = (Aa)xb = ax(Ab) 

and so we can write these equivalent expressions in the form Aa x b. 
At this stage we wish to establish the non-trivial relation 

ax(b + c) = axb + axc. 

As a first step let (3 (Fig. 15) be the projected vector of b on a plane 
p perpendicular to a. We have j8 = b sin 6 and so 

axb = ax (3. 
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Now let y be the projected vector of c on the plane p. Then |3 + y is 
the projected vector of b + c on the plane p. Hence 


and 


axe = ax y 
ax(b + c) = ax(p+y). 



Fig. 15 


Since a is orthogonal to (3, it follows that a x (3 is a vector in the 
plane p with magnitude af3 in a direction orthogonal to (3. That is, 
a x |3 is obtained by rotating the vector a$ through a right-angle in the 
plane p. Similarly a x y is the vector obtained by rotating the vector 
ay in the plane p through a right-angle in the same sense as the 
previous rotation. Thus ax(3 + axyis the sum of the two vectors aft 
and ay followed by the appropriate rotation in the plane p through a 
right-angle. This sum «(P+y) rotated through a right-angle in the 
plane p is the vector a x (|3 + y) and so 

ax(P + y) = axp + axy. 

Accordingly, we have 

ax(b + c) = axb + axc. 

As a consequence we may remove brackets as in ordinary algebraic 
multiplication. That is, vector multiplication is distributive with 
respect to addition. 

Let the vectors a and b have components a ti a jt a k and b ir b jy b k 
with respect to a fundamental system i, j and k. Then 

axb = (fl,i+ cj + a k k) x (b t i+b J + b k k). 
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Removal of brackets and use of the vector product properties (9.1) of 
the vectors of the fundamental system yield 

axb = (a j b k -a k b J )i + (aJ> i -a i b k )j + (a i b J -a J b i )k. (9.2) 

This expression may be written in the determinantal form 

flj &} Afc 

b, b, b k . 

i j k 

Example. Let A , B and C be given by the vectors a, b and c. Show 
that ^(bxc + cxa + axb)isa vector orthogonal to the triangle ABC 
with magnitude equal to the area of the triangle. 

The vectors AB and AC are given by b — a and c — a. Hence the 
required vector is ^(b — a) x (c — a) from which the required result 
follows. 

Exercises 

1. Prove that (axb) 2 = a 2 b 2 — (a.b) 2 . 

2. Given a = 2i—3j + k, b=—i + k and c = 2j —k, (i) calculate axb, 
b x c and c x a in terms of i, j and k, (ii) obtain a unit vector orthogonal to 
both b and c, (iiij"calculate the area of the parallelogram with the diagonals 
a+bandb + c. 

3. Given a.b = a.c, axb=axc and a is not the zero vector, show 
that b = c. 

Let a, b and c represent the vectors BC, CA and AB respectively. 
Show that axb=bxc = cxa and deduce the sine rule for a triangle. 

10 Scalar Triple Product 

We now consider a. b x c. There is no ambiguity in this expression 
as the vector product b x c must be evaluated and the scalar product of 
this vector then taken with a. The result is called the scalar triple 
product of the three vectors. 



Fig. 16 
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The vector b x c is the vector be sin 6 n where 6 is the angle 
(0 < 6 n) between b and c, and n is the unit vector orthogonal to both 
b and c such that b, c, and n form a right-handed system (Fig. 16). 
The magnitude be sin 6 of the vector b x c can be interpreted 
geometrically as the area of the parallelogram formed by b and c. 
Let a be the angle (0 < a ^ 77 -) between the vectors a and n. Then 
a .b x c = abc sin 6 cos a. 

The volume of the parallelepiped formed by the three vectors a, b 
and c is equal to the product of the area of the parallelogram formed 
by b and c with the height of the parallelepiped. This height is the 
component of a in the direction n and so is a cos a. Accordingly, the 
volume of the parallelepiped is abc sin 6 cos a. 

That is, the scalar triple product a. b x c equals the volume of the 
parallelepiped formed by the three vectors when set off from a 
common point. 

If a.bx c is positive, cos a is positive and so a. is acute. Hence n 
and a extend on the same side of the plane formed by b and c. How¬ 
ever, if a.bx c is negative, a lies between 77-/2 and 7 r and so n and a 
extend on opposite sides of this plane. 

We extend the definition of right- and left-handedness to non- 
orthogonal sets of three vectors thus: the vectors a, b and c are 
right-handed (left-handed) if a and n are on the same side (opposite 
sides) of the plane formed by b and c. The reader is asked to verify 
that when a, b and c are mutually orthogonal, this definition coincides 
with that of Section 7. 

That is, a, b and c form a right-handed or a left-handed system 
according as a.bxc is positive or negative respectively. 

By a similar argument the volume of the parallelepiped is given 
by b.cxa and c.axb and these triple products are positive or 
negative according as a, b and c form a right-handed or left-handed 
system respectively. Hence we have 

a.bxc = b.cxa = c.axb. 

Since c.axb = axb.c, it follows that 

a.bx c = ax b.c. 

That is, we may interchange the dot and cross in a scalar triple 
product. 

Further, we see from the left-handed or right-handed properties 
that 

a.bxc = —b.axc. 

The reader is asked to verify that in all cases the scalar triple product 
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of a, b and c taken in the cyclic order a, b, c equals a.b x c but equals 
— a. b x c for any anti-cyclic order of the vectors. 

From the geometrical interpretation it follows that the scalar 
triple product a.b x c vanishes if one of the vectors is the zero vector 
or the three vectors a, b and c are coplanar. 

Conversely, let the three non-vanishing vectors a, b and c satisfy 
the relation a.bxc = 0. Then a is orthogonal to bxc. But by the 
definition of a vector product, both b and c are orthogonal to b x c 
and so a, b and c are all orthogonal to b x c. Accordingly a, b and c 
are coplanar vectors. 

Let a, b and c have components a u a Jt a k \ b u b jy b k and c u c jt c h 
respectively with reference to the fundamental system i, j and k. 
Then by (9.2) we have 


b x c = (bjC k - Vi)* + (Vi - Vfc)j + (V “ Vi) k 

and so 

a . b x c = a$f k - V/) + Vi - V) + W ~ V)- 
This result may be written in the determinantal form 


Q-j Cl k 


a.bxc = 


bi bj b k . 

c i c j c k 


We agree to use the notation 

[abc] = a.bxc. 

We note that the vectors of a fundamental system satisfy the 
relations 


[ijk] = [jki] = [kij] = 1, 
[ikj] = [kji] = [jik] = -1. 


Exercises 

v 1. Obtain the value of A which makes the vectors i—j + k, 2i+j — k and 
Ai — j + Ak coplanar. 

2. If a, b and c are orthogonal vectors, show that 

[abc] 2 = a 2 b 2 c 2 . 

3. Show that 

(a+b).(b+c)x(c+a) = 2[abc]. 

11 Vector Triple Product 

Consider the vector triple product ax (bxc). This vector is 
orthogonal to both a and bxc. But b x c is orthogonal to both b and 
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c. It follows that the vector a x (b x c) lies in the plane of b and c and 
so we may write 

ax(bxc) = Ab+/zc. 

Since a x (b x c) is orthogonal to a, its scalar product with a is 
zero and so 

Aa.b+/xa.c = 0. 


Accordingly we may write 

A = pa.c, fi = —/>a.b, 

where p is still undetermined. Thus we have 

ax(bxc) = />[a.cb — a.be]. (1) 

First we consider the case when a = b and take the scalar product 
of both sides of this equation with c to obtain 

c.bx(bxc) = p[(b.c) 2 —b 2 c 2 ]. 

Interchanging the dot and the cross of the scalar triple product we 
have 

— (bxc) 2 = p[(b. c) 2 —b 2 c 2 ]. 

Let the angle between b and c be rp. Then this equation yields 

— b 2 c 2 sin 2 tp = p[b 2 c 2 cos 2 (p—b 2 c 2 ], 

from which p= 1 in this special case. That is, 

bx(bxc) = b.cb —b 2 c. (2) 

Now return to the general case when a and b differ and take the 
scalar product of both sides of equation (1) with b to obtain 

b.ax(bxc) = p[a.cb 2 —a.bb.c]. 

Interchange of the cyclic order of b, a and bxc followed by substitu¬ 
tion from equation (2) yields 

— a.[b.cb — b 2 c] = />[ a.cb 2 —a.bb.c], 

from which we see that p is also unity in the general case. That is, 

ax (bx c) = a.cb — a.be. (11.1) 

Note carefully that (a x b) x c is not in general the same vector as 
ax (bxc) because (axb)xc is a vector in the plane of a and b. 
Accordingly it is essential to retain brackets in a vector triple product. 
We have 

(a x b) x --cx(axb) 

= a.cb — b.ca. 
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Exercises 

1. Establish the formula for a x (b x c) by introducing the components 
of the vectors with respect to a fundamental system of vectors. 

, 2. Show that i x (j x k) = 0. 

^ 3. If a=i + 2j — 3k, b=2i—j + kand c= — i + 4j—k, calculate (i) a.bx c, 
(ii)(a x b) x c, (iii) b x (a x c), (iv) (a x b) x (c x a), ( v )[a x b b x c c x a], 
(vi) (a x b) x {(bxc)x(cxa)}. 

v 4. Show that a x (b x c) = (a x b) x c if and only if b=0 or c is parallel 
to a or b is orthogonal#to both c and a. 

,5. If a is not the zero vector, prove that every solution of the equation 
a x x=b is given by x=(b x a)/a 2 + Aa, where A is an arbitrary scalar. 

12 Products of Four Vectors 

There is no ambiguity in the expression axb.cxd as the two 
vector products must be formed before the scalar product can be 
evaluated. The interchange of dot and cross yields 

axb.cxd = a.bx(cxd) 

= a.(b.dc — b.cd) 

= a.cb.d—a.db.c. (12.1) 

In particular, choose c = a and d=b to obtain 

(a x b) 2 = a 2 b a —(a.b) 2 . (12.2) 

Next let us consider (a x b) x (c x d) and evaluate as the vector 
triple product a x (c x d), where a = a x b. The result is 

(axb)x(cxd) = [abd]c — [abc]d. 

The calculation can also be made from (a x b) x (3, where (3 = c x d 
with the result 

(axb)x(cxd) = [cda]b — [cdb]a. 

Equate the two expressions for (a x b) x (c x d); then solve for d, 
provided that [abc] is not zero, to obtain 

A [bcd]a + [cad]b + [abd]c , 10 ^ 

a ~ [abc] 

This relation shows how any vector d can be expressed as a linear 
combination of any three given non-coplanar vectors a, b and c. 
That is, any four vectors are always linearly dependent. 

Exercises 
1. Prove 

(i) ax(bxc)+bx(cxa) + cx(axb) = 0 




1 ] 


RECIPROCAL BASIS 


21 


(ii) [bxc cxa axb] = [abc] 2 , 

v(iii) (b x c) x (a x d) + (c x a) x (b x d)+(a x b) x (c x d)= — 2[abc]d, 

(iv) b xc.ax d+ cxa.bx d+ axb.cxd=0, 

(v) b x ca. d+c x ab. d+a x be. d= [abc]d. 

2. The unit vectors a, b and c define a spherical triangle ABC on a unit 
sphere with centre at the origin O. Let a, (3 and y denote the arcs of the 
great circles forming the triangle, labelled so that a is opposite to the 
angle at A, etc. Show that (a x b) x (a x c) = [abeja and deduce that 
sin a sin (3 sin C=[abc]. Hence obtain 

sin A _ sin B _ sin C 
sin a sin ft sin y 

This relations is known as the sine law for spherical triangles. 


13 Reciprocal Basis 

A set of three non-coplanar vectors e lf e 2 and e 3 is called a basis 
because any vector can be expressed linearly in terms of any three 
non-coplanar vectors by equation (12.3). The basis is said to be right- 
handed or left-handed according as [e 1 e 2 e 3 ] is positive or negative. 

A second basis e 1 , e 2 and e 3 f is said to be reciprocal to e l5 e 2 and 
eo if 


(13.1) 


The vector e 1 is orthogonal to both e 2 and e 3 and so we may write 
e 1 = Ae 2 x e 3 . Substitution in e 1 .e 1 = 1 yields A^e^g] = 1. Hence we 
have 


ei.e 1 = 1, 

ei.e 2 

= 0 , 

ej.e 3 = 0, 

ea.e 1 = 0, 

e 2 .e 2 

= 1, 

e 2 .e 3 = 0, 

esj.e 1 = 0, 

e 3 .e 2 

= o, 

e 3 .e 3 = 1. 


»i _ 


Go x ec 


e 2 = 


eq x e. 


e° = 


e! x e. 


[ ei e 2 e 3 ]' ~ [e^eg]' ~ [e ie2 e 3 ] 

From the relative symmetry in e 1} e 2 , e 3 and e 1 , e 2 , e 3 we have 
similarly that 


e, = 


e'xe° 

[e 1 e 2 e 3 ]’ 


e«> = 


e° x e A 
[e 1 e 2 e 3 ]’ 


e s = 


e A x 

[e 1 e 2 e 3 ] 


(13.3) 


Exercises 

•' 1. Obtain a set of vectors reciprocal to the three vectors — i+j + k, 
i— j + k and i+j—k. 

f The indices 1, 2 and 3 of e 1 , e 2 and e 3 merely serve as distinguishing 
‘labels’ and do not possess any significance as power indices. 
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v 2. If e lf e 2 e 3 form a self-reciprocal basis, show that they are the unit 
vectors of a fundamental system. 

3. Prove that (i) [e 1 e 2 e 3 ][e 1 e 2 e 8 ] = l, 

(ii) e 1 xe 1 + e 2 xe 2 +e 3 xe 3 = 0, 

(iii) [e 1 e 2 e 3 ] 3 [e 1 e 2 e 3 ] = [e 2 x e 3 e 3 x e x e x x ej. 

Miscellaneous Exercises 

1. If b x c is not the zero vector, obtain x in terms of a, b and c given 
that a.x=l, b.x=0 and c.x=0. 

2. If a, b and c are three non-coplanar vectors, show that the vectors 
b x c, c x a and a x b are non-coplanar and prove that any vector d can be 
expressed in the form 

d = {a.dbxc+b.dcxa+c.dax b}/[abc]. 

3. Show that the internal bisector of an angle of a triangle divides the 
opposite side in the ratio of the other two sides. Deduce that the internal 
bisectors of the angles of a triangle are concurrent. 

4. The coplanar straight lines AL, BM and CN are concurrent. Prove 
that the points of intersection of AB and LM, BC and MN, CA and NL 
are collinear. 

5. Show that the mid-points of the diagonals of a complete quadrilateral 
are collinear. 

6. Show that the sum of the four vectors of magnitudes equal to the 
areas of the faces of a tetrahedron and directed outwards perpendicular to 
the faces is zero. 

Generalize this result to a closed polyhedron. 

7. If the vector triple e 1 , e 2 , e 3 is reciprocal to the triple e x , e 2 , e 3 , 
show that any vector a satisfies the relation 

a = (a.e 1 )e 1 +(a.e 2 )e 2 +(a.e 3 )e 3 

and 

a = (a.e 1 )e 1 +(a.e 2 )e 2 +(a.e 3 )e 3 . 

8. Show that (axb)*(axc)=(b.c)(a.a)-(a.c)(a.b). Hence, with the 
notation of exercise 2, page 21, deduce the cosine law for spherical 
triangles 

cos a = cos j8 cos y+sin j8 sin y cos A. 

9. Obtain the angle (a) between two faces, (b) between a face and an edge 
which intersects it, of a regular tetrahedron. 

10. D , E and F are points on the sides BC, CA and AB of a triangle 
ABC such that AD, BE and CF are concurrent. Use vector methods to 
prove Ceva’s theorem that 

BDCEAF 

dc'ea'fb ~ L 
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Applications to Space 

Geometry t 

In this chapter vectors are used in two different ways. Firstly, the 
position of a point in space is denoted by a vector whose initial point 
is fixed. The other end of the vector determines the point uniquely in 
space. Such vectors will be denoted by bold Roman letters. Secondly, 
a vector may merely determine a direction in space irrespective of its 
initial point. Such vectors will be denoted by bold Greek letters. 

14 Straight Line 

There is a unique straight line (Fig. 17) through A, given by the 
vector a referred to some fixed initial point O, called the origin, in 
the direction of the vector a. Let r be the position vector of any point 



o 


Fig. 17 

R on this straight line. From OR = OA+AR we have 

r = a + tta 

f The content of this chapter is not required for the understanding of the 
remainder of the text. 


3—v.a. 
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where the value of the parameter u depends on the position of R on 
the straight line. As u varies from —oo to +oo the point R traces out 
every point of the straight line. 

We have AR = r — a and so AR 2 = (r—a) 2 = u 2 a 2 . Accordingly, u 
measures the actual distance AR if and only if a is chosen to be a unit 
vector. 

The parameter u can be eliminated from the equation r — a + ua by 
writing it in the equivalent form 

(r —a)xa = 0. 

Next we calculate the perpendicular distance (Fig. 18) BM from 
B t given by the vector b, to the straight line r=a + ua. We have 

AB=b—a and \AB x a| — aAB sin 6 — aBMy where 8 is the angle 
BAM. 



Fig. 18 


Accordingly, 

BM 2 = {(a—b) xa} 2 /a 2 . 

An alternative method yields 

BM 2 = AB 2 — AM 2 

= (a—b) 2 —{(a — b). a/a} 2 
_ (a—b) 2 a 2 —{(a —b).a} 2 


The equality of the two expressions for BM 2 is verified from equation 

( 12 - 2 ). 

Example. Show that if the two straight lines r=a+ua and r = h + v$ 
intersect then (a—b).ax(3 = 0 but a x(3#0. 
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If the straight lines intersect, values of u and v exist for which 
a + «a=b + ©p. The required result is obtained by taking the scalar 
product with a x (3. Note that the vanishing of a x (3 implies that the 
straight lines are parallel. 

15 Plane 

Consider the plane (Fig. 19) through A, given by the vector a, 
containing the directions given by the vectors a and (3. By the 
parallelogram law of addition of vectors, the vector AR corresponding 



to any point R on the plane of a and (3 is given by wa+^P where u 
and v are parameters. Then the position vector r of R is given by 

r = a + wot+ ©(3. 

This equation represents the required plane. To each point R on it 
there corresponds a pair of values of w and v. 

We may eliminate the parameters u and v to obtain the equation 
of the plane in the form 

(r—a).otx(3 = 0. 

This equation can be obtained directly by noting that ax (3 is a 
vector orthogonal to the plane and r—a is a vector lying in the plane. 

The equation of the plane through the three points A, B and C, 
given respectively by the vectors a, b and c, is obtained by noting 
that the vectors b — a and c — a are parallel to the plane. Hence the 
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equation satisfied by r, the position vector of any point R on the plane, 
is 


That is, 
where 


r = a + w(b — a) + ©(c — a) 
= (1 — u — ©)a+wb + vc. 

r = Aa+pb + vc, 


X + fJi + V = 1 . 


Next we calculate the perpendicular distance BL (Fig. 20) from B , 
given by the vector b, to the plane (r-a).(oxp) = 0. We note that 



BL is the component of the vector BA = a—b in the direction of BL. 
This direction is parallel to the vector a x p. Hence we have 

_ ±(a-b).axp 
|axp| 


16 Shortest Distance between Two Skew Lines 

Consider (Fig. 21) the two straight lines through A and B , given by 
the vectors a and b, in the respective directions a and p. The corre¬ 
sponding equations are r=a + ua and r=b+©p. The shortest distance 
between a point and a straight line is the perpendicular to it from the 
point. Consequently the shortest distance between two skew lines is 
along the line which is perpendicular to both of them. That is, the 
shortest distance between these two straight lines is LM, where LM is 
perpendicular to both the vectors a and p. Accordingly, the vector 
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a x p is along LM. But LM is the component of the vector AB in the 
direction of LM and so we have 


LM = (a-b^oxp 

l«xPI 


since (a x P)/|a x (3| is a unit vector in the direction of LM. 



o 


Fig. 21 


We deduce that the necessary and sufficient condition that the two 
straight lines intersect is (a—b). a x (3 = 0 provided that a x (3 ^ 0. We 
note that this condition agrees with the result in Section 14. 


Miscellaneous Exercises 

1 . Calculate the perpendicular distance from the point (—2, 1, 5) to the 
straight line joining the points (1, 2, —5) and (7, 5, —9). 

2 . Calculate the perpendicular distance from the point ( — 2, 3, —1) to 
the plane determined by the three points ( — 4, 4, —8), ( — 5, 8, —11) and 
(-7,-2,-8 ). 

3. Calculate the shortest distance between the straight lines joining the 
points (—14, 8, 6), ( — 11, 4, 1) and the points (3, 5, 5), (6,11, 8). 

4. Show that the straight line through the point represented by the 
vector b perpendicular to and intersecting the straight line r=a+«a is 
given by r=b+mx{(a-b) xa}. 

5. Show that the equation of the plane through the points A, B and C 
given by the vectors a, b and c respectively can be written in the symmetrical 
form 


(3r—a—b—c).(bxc+cxa+axb) = 0. 
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6. Show that the point of intersection of the straight line (r — a)xa=0 
and the plane (r — b) .(3 = 0 is at the point given by the vector a+Aa, where 

v _ (*>-»)•P 

a.p 

7. Show that the line of intersection of the two planes (r — a).a=0 and 
(r — b) .p = 0 is given by (r—c) x (a x (3) = 0, where 

c(axp) 2 = a.a[px(axp)]+b.p[ax(pxa)]. 

8. The Plucker vectors of the straight line through the point A given by 
the vector a and parallel to the direction of the vector a are p and q, where 
p = a/|a| and q=axp. Show that 

(i) p is independent of the choice of the point A on the straight line, 

(ii) the perpendicular distance from the point given by the vector b from 
the straight line is |pxb + q|, 

(iii) the straight lines corresponding to the Plucker vectors p 1? q x and 
p 2 , q 2 intersect if and only if pi.q 2 + p 2 .qi = 0 and PiXp 2 #0, 

(iv) the equation of the plane through A and the line with Plucker vectors 
Pi, qi is 

r.qi-[rapj = a.q^ 
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Differential Vector Calculus 

17 Derivative of a Vector 

Let p(w) denote a vector function of the scalar variable u, where u is 
restricted to some range of values u x < u < u 2 . By this statement we 
mean that p(w) is uniquely determined when u is given a value in its 
range. 

As an example consider the vector p = a + wot. As u varies, p 
represents the vectors joining the origin to any point of the straight 
line through the pointf a in the direction a. For this vector, u may 
range from — oo to +oo. 

The vector function p (m) is said to be continuous at u 0 if, given 
any positive number e, it is possible to find another positive number 
S such that |p(«) — p(«o)| <c when \u — u 0 \ < S. This condition is also 
written in the form lim p(«)=p(«o)- 

The derivative p(w) is defined by 

p(«) = 4 = Um P(«+^)-P(«) 

du Au-* 0 AU 

provided that this limit exists. In this case we also say that p(w) is 
differentiable. (The dot above the letter will always denote 
differentiation with respect to u.) 

Since p (u) is itself a vector function of w, we can consider its 
derivative with respect to u. If this derivative exists we denote it by 
p(w), where 

m - ^m) - 3 - 

Similarly, higher derivatives of any order can be defined. 

f When the vector a represents the directed line segment OA, where O is 
the origin, we occasionally refer to A as the point a. 
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Example. Obtain the derivative of the constant vector c. Let p(w) = c 
and so p(u+Au) = c from which £(w) = 0. That is, the derivative of a 
constant vector is the zero vector. 


18 Derivative of a Sum of Vectors 
Divide both sides of the identity 

{p(tt +Au) + q(w +Au)}— (p(m) + q(w)} 

= {p(w +Au)— p(m)} + {q(w +Au) — q(w)} 

by Au and proceed to the limit as Au tends to zero. It follows that 



It follows from the commutative and associative laws of addition 
of vectors that for a finite number n of vectors p,(tt) we have 


Hi-)-it 


19 Derivative of the Product of a Scalar and a Vector Function 

Consider p(u) =f(u)a(u), where f(u) is a scalar function of u and 
a( m) is a vector function of u. Then we have 


p {u+Au) — f(u+Au)a(u+Au). 

Divide both sides of the identity 
p(u+Au)—p(u) 

= f(u +Au){a{u+Au) — a(w)} + { f(u+Au) —f(u)}a{u) 
by Au and proceed to the limit as Au tends to zero. It follows that 

d ,. . r da df 


Exercises 

1. If n, a and b are constant and p = a cos nt+h sin nt, prove that 
px^ = rcaxb and ~^+n 2 p = 0. 


dt 2 


2. Find p satisfying the equation p(«) = a«+b, where a and b are 
constant vectors, given that p and p both vanish when «=0. 

3. Show that ^{vlp)=vlp-pvlp 2 - 
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20 Derivative of a Scalar Product 

Consider / (u) — a(w). b(w), where a and b are vector functions of 
the scalar u. Then we have f{u+Au) — a{u+Au).h{u-\-Au). Divide 
both sides of the identity 

f(u+Au)~f(u) 

= a (u +Au) . {b(« +Au) — b(w)}+ {a(u+Au) — a(«)}. b («) 
by Au and proceed to the limit as Au tends to zero. It follows that 



db da. , 
a ’ du +du' 


In particular, we have 




But a 2 = a 2 and so 



2 a —• Hence we have 
du 


da 


a. 


du 



Further, if a is a vector of constant magnitude , = 0 and so a. =0. 

du du 

That is, the derivative of a vector of constant magnitude is either 
orthogonal to it or zero. 


21 Derivative of a Vector Product 

Consider p(w) = a(w) x b(w), where a and b are vector functions of 
u. Then \*(u+Au) = a(u+Aii) x h(u-\-Au). Divide both sides of the 
identity 

p(tt+dtt)—p(tt) 

= a (u +Au ) x {b(«+d«) — b(t/)}+{a (u +Au ) — a(w)} x b(w) 
by Au and proceed to the limit as Au tends to zero. It follows that 


d . , . dh da , 

-j-faxb) = ax-r-i-j-xb. 
du ' du du 


It is important in this case not to interfere with the order of the 
vectors as vector multiplication is anti-commutative. 
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We immediately deduce that 

-^-(a.bxc) = a.bxc+a.^-(bxc) 

du v ' du 

= S.bxc + a.bxc+a.bxc. 

Similarly we obtain 

^-{ax(bxc)} = ax(bxc) + ax(bxc)+ax(bxc). 


Exercises 

1. Show that (i) ^ (pxp) = pxp, (ii) ^ [ppp] = [PPP]• 

2. If a is a unit vector in the direction of b, show that a x a = (b x b)/(b.b). 

3. If ^=uxa and ^wxb, show that 

j du du 

^(axb) = c*) x (a x b). 

dx • • 

\/ 4. If r x —= 0, show that r has a fixed direction. (Hint: Let x—f(u)e(u), 
du 

where e(w) is a unit vector. Deduce that e x e = 0 and e.e = 0.) 


22 Taylor’s Theorem for a Vector Function 

Suppose that the vector p (u) and its derivatives up to order n— 1 
are continuous for a ^ u ^ a + h and p (n) (w) exists for a<u<a+h. Under 
these conditions we prove Taylor’s theorem that 

p (a + h) = p(tf) + hp\a) H- 

+ 7 5-rr-r ll n ~ 1 p (n " 1) (fl) + ”7 h n p in \a + 6fl), (22.1) 

(n— 1)! n\ 

where 0< 6< 1. 

If, in addition, lim {(h n ln\)p in \a+6h)} = 0 the vector function 

71-» CO 

p (a+h) has the Taylor expansion 

p (a + h) = p(a) + hp'{a) H- h n p in) (a) H-. 

Define the vector q(w) for 0 ^ u ^ h by 
q(w) = p {a + u) - p(«) - up'{a) -^ 
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where the constant vector X is chosen so that q(/z) = 0. We verify that 
q(0), q'(0), q"(0), . •q (n-1 ) (0) 

are all zero. 

Now choose an arbitrary vector c and define the scalar Q(u) by 

Q(u) = c. q(w). 

It follows that 

Q( 0), QX 0), 0"(O), • - ^-^(O) and Q(h) 

are all zero. 

Successive application of Rollers theorem, yields 

(1) Q\h^) = 0, where 0 <h x <h since <2(0) = 0 and Q{h) = 0, 

(2) Q"{h 2 ) = 0, where 0 <h 2 <h x since <2'(0) = 0 and Q\h i) = 0, 


(«) Q {n \h n ) = 0, where 0 <h n <h n _ 1 since ^ (n “ 1) (0) = 0 and 

Q^-'XK- i)=o. 

Since 0 < h n < h n _ 1 < • • • < h x < h, there exists a number 6 for which 
h n = 6h and 0 < 6 < 1. But 

Q n (u) = c.q (n) («) = c. [p (n) (« + w) — X] 

and so 

c[p< w >(«+0/z)-X] = 0. 

However, c is an arbitrary vector and so X = p (n) (a + 6h). Substitu¬ 
tion in q(/z) = 0 yields equation (22.1) as required. 

23 Derivative of a Vector Referred to a Fundamental System 

Let the vector p(w) have components p { {u), p/u) and p k (u) referred 
to a fundamental system i, j and k. That is, 

p(w) = pt(u)i +pj(u)j +p k (u)k. 

Then we have 

W-SW) + 5M+S(^ 

Since i, j and k are constant vectors, we have 

P(«) = Pii+Prf+PjM. 


Exercises 

1. Given that r =a cos u i +a sin u j + buk, show that r 2 — a 2 + b 2 , 
(r x r) 2 = a 2 (a 2 +b 2 ) and [r if *r] = a 2 b. 
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2. If a=«i—« 2 j+« 3 k and b=sin u i+cos «j, calculate (i) ^ (a. b), 


24 Partial Derivatives of Vectors 

Let p (u, v) denote a vector function of the independent scalar 
variables u and v. By this statement we mean that p(w, v) is uniquely 
determined when u and v are given values within their ranges. 

As an example consider the vector p = a + ua + ©(5. As u and v vary 
p represents the vectors joining the origin to any point of the plane 
through the point a in the directions given by a and (3. 

The partial derivative p u is defined by 


Vu = 


dp 

8u 


Hm p (u+Au y v)-p{u , v) ^ 

All-*0 


provided that this limit exists. 

Similarly, the partial derivative p„ is defined by 

__ 1im p(u,v + Av)-p(u,v) 

Vv ~ Sv ~ £?. Av 


provided that this limit exists. 

More generally, p may denote a vector function of any number of 
independent scalar variables and it is clear how partial derivatives 
are defined. 

We may in the usual way define partial derivatives of higher order. 
It may be shown, as in real analysis, that under suitable conditions 
partial differentiation is commutative. For example, we have 

8 dp _ 8 dp 
8v 8u 8u 8v 


As in the preceding sections we can similarly establish the results 
that 


8 _ 

8u 


(p+q) 


_ 8p 8q 
8u + 8u’ 


d_ 

du 


-da df 


a, 
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8 ,. 8b 8a , 

^(a.b) = a.^+- .b. 


8 . .. 0b 2a , 

^(axb^ax^+^xb. 


+ ^ + ^> = t i+ ^ + t k - 


Note carefully that the indices j, / and k refer to the components 
of a vector with respect to a fundamental system. They will not be 
used as partial derivative indices. 


Exercises 

,1. If p = e m ' i+(2u—v)]+v sin u k, calculate (i) (ii) (iii) ^ 

' du ' ' du 2 v ' dudv 

r \ 8^ f \8p dp 


2 . 


If a=i«wi+M«> 2 j—w 3 k and b=« 3 i—«v«j+M 2 wk, calculate 


(0 


g 2 a 

dv 


d 2 a d 2 b 

at the origin, (ii) ^ x ^ at the point (1, 1, 0). 


CHAPTER 4 


Applications to Differential 

Geometry t 

25 Curve and Tangent Vector 

Consider the vector function r(w) of the scalar variable u. To each 
value of u there corresponds a vector whose initial point is at the 
origin. The other end of the vector determines a point in space. As u 
varies, this point traces out a curve. To avoid having the curve 
degenerate to a point, the case when r(w) is a constant vector is 
excluded. 



f Sections 26 and 27 may be omitted by the reader who is not interested in 
differential geometry. The remaining sections are vital to the subsequent 
development. 


36 
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Let P and Q (Fig. 22) correspond to values u and u+Au of the 

scalar variable. Then OP and OQ are given by the vectors r(w) and 
x(u+Au) respectively. 

The tangent at a point P on the curve is defined to be the limiting 
position of the chord PQ as Q tends to coincidence with P. We have 

PQ = OQ — OP — x(u+Au) — r(u) 

and so 

PQ _ t(u+Au) — t(u) 

Au Au 


Let Q tend to P. Then Au tends to zero and we obtain 


lim = r(u). 
^u-*o Au 


PQ . 

But lim is a vector in the direction of the tangent and so the 
au-*o Au 


vector r(u) is a tangent vector to the curve r(u) at the point corre¬ 
sponding to the scalar variable u. 

The vector r(w) is not tied to a fixed point as an initial point and 
so by the convention introduced in Chapter 2 we ought to represent 
it by a Greek letter. However, in this and in the following chapters we 
shall not employ this convention. It is merely necessary to remember 
that r( m) is tied to the origin O but other vectors are free in space. 

Let the coordinates of the point P given by r(w) be x } y and z 
referred to a set of rectangular axes with O as origin. Further, let 
i, j and k denote the fundamental system of vectors with respect to 
these axes. Then we have 


r (u) = x(u)i +y(«)j + #(tt)k, 

where x(u), y(u) and z(u) are scalar functions of u. 

Introduce the scalar s which measures the distance along the curve 
from some fixed point of the curve. Then we may select s in place of 
the parameter u to obtain 

r (s) = jc(s)i +y(s)j + #(s)k. 

The vector r'(s), where the prime sign denotes differentiation with 
respect to s , is now a tangent vector. 

From ds 2 = dx 2 + dy 2 + dz 2 i we have 

W(s)] 2 + [y'(s)] 2 + [z'(s)] 2 = 1. 
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Accordingly 

r \s ). r'(s) = [*' (*)i +/(*) j + z'(s)k] 2 

= W(s)] 2 +[y'(s)] 2 +[A*)] 2 
= 1 . 

That is, 

T( S ) S r' (s ) = * 

is the unit tangent vector to the space curve r(s) and it clearly points 
in the direction of increasing s. 

Let A and B be points on the curve r(u) corresponding to the 
parameter values u A and u B . If r(w) is continuous for u A ^u^u B , we 
say that the curve joining AB is a regular arc. A double point on a 
curve is a point at which the curve intersects itself or is tangent to 
itself. Further, a regular curve consists of a finite number of regular 
arcs without double points. That is, there can at most be a finite 
number of points with discontinuous tangents on a regular curve. 

26 Frenet Formulae 

The vector T(s) = r'($) is the unit tangent vector to the curve r($), 
where s measures arc-distance along the curve from some fixed point 
on it. Differentiation of r'(s). r'(s) = 1 yields r'(s). r"(s) = 0. That is, 
the vector r"(s) is orthogonal to r'(s). The direction determined by 
r"(s) is called the positive direction of the principal normal. Let 
N(s) denote the unit vector in the positive direction of the principal 
normal. Accordingly a function #c, positive by convention, called 
the curvature exists such that 

N(*) = ir”( S ), (k>0). 

This equation can be written in the form 


The unit binormal vector B(s) is defined to be the unit vector 
orthogonal to both the tangent vector T and the principal normal 
vector N such that T, N and B form a right-handed system of unit 
vectors. Accordingly 

B = TxN, N = BxT, T = NxB. 

B is a unit vector and so B' is orthogonal to B (see Section 20). 
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4] 

Further we have B.T=0 and so B'.T+B.T' = 0. But B.T' = B.icN=0 
and so B'.T=0. Hence B' is also orthogonal to T. Accordingly a 
function r($), called the torsion, exists such that 


B' = 


dB 

ds 


-tN. 


It is merely a matter of convention that r is introduced by this equation 
with a negative sign. Note carefully that the sign of r itself cannot be 
fixed by convention to be positive, but is positive if B' and N have 
opposite directions. 

We have 

N = BxT 

and so 

N' = B'xT+BxT' 

= — tN x T+B x kN 


= rB — kT. 


The three equations for the derivatives T', N' and B' in terms of the 
basis formed by T, N and B are called the Frenet formulae. Col¬ 
lecting the Frenet formulae together, we have 


T' = kN 

N' = -kT +tB ► 

B' = -rN. 


(26.1) 


Introduce the Darboux vector co = rT+ kB. Then the Frenet formula 
can be exhibited in the symmetrical form 


T' = to x T, N' = to x N, B' = co x B. 


The three vectors T, N and B are called the trihedral at a point. 

Example 1. Show that the necessary and sufficient condition that a 
curve be a straight line is that k = 0 . 

If the curve is a straight line, the vector T is constant and so 
T' = 0. Hence k = 0. 

Conversely, given k — 0, we have T'=0 and so T is a constant vector 
and the curve is a straight line. 

Example 2. The necessary and sufficient condition that a curve lie in 
a plane is that r=0. 

If the curve lies in a plane, select the origin in this plane. Both r' 
and r 9 are vectors lying in this plane. Hence T and N lie in the plane 
and so B is a constant vector, its direction being orthogonal to the 
plane. Accordingly B' = 0 and this implies that r = 0. 
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Conversely, given r = 0 we have B' = 0 and so B is a constant vector. 
Let the curve be given by r($). Then we have 

^(r.B) = T.B + r.B' = 0 
as v 7 

and so r.B is a constant. That is, the vectors from the origin to points 
of the curve are all orthogonal to a fixed direction B. Hence the curve 
is a plane curve. 


27 Curvature and Torsion 

We have 

— T, 

r" = T' = #cN, 

and so 


x m = fc'N + fcN' 


Thus 

and 


= #c'N + #c( — #cT + rB) in virtue of (26.1) 
= -fc 2 T+fc'N + fcrB. 

r". r" = K 2 

r'.r'xr" = T.kNxktB = k 2 t. 


Hence we have 

k = and r 


r .r xr 


r .r 


(27.1) 


In general the parameter u in r(w) is not the arc-distance s and so 
we cannot evaluate k and r directly from these equations. 

For any vector a we have 


. _ da 
a du 


da ds 
ds du 


Using this result in conjunction with the Frenet formulae (26.1) we 
have 

r=iT, T = iicN, N = s( - kT + rB), B = - irN, 
r = fT-|-f 2 KN, 

r = s T + sskN + 2sskN + s 2 k N+ s 3 k( — kT + rB) 

and so 

r = (s — s 3 k 2 )T + (3ssk + i 2 /c)N + 5 3 /ctB. 
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Thus we obtain 


r x r = s 3 kB, 

|rxr| = |$ 3 [/c = [r| 3 «r, 
r.rxr = s 6 k 2 t. 


rxr 


and r 


r. if x r 
|rxr| 2 


Example 1. Calculate the curvature and torsion of the curve 
r = a(3u — w 3 )i + 3aw 2 j + a(3w + w 3 )k. 

We have 

x = 3«(1 — M 2 )i + 6auj + 3«(1 + « 2 )k, 

? = — 6awi + 6aj + 6flwk, 
r = — 6oi +6<zk, 

Accordingly, calculation yields 

|f| 2 = 18a 2 (l +M 2 ) 2 , 

rxr = 18 o 2 (m 2 —l)i —36a 2 wj +18 o 2 (m 2 +l)k, 
|rxr| 2 = 2.18 2 o 4 (1 + m 2 ) 2 , 
r.ifxr = rxr.r = 216a 3 . 


(27.2) 


Thus we obtain on substitution into equations (27.2) that 

1 

* T 3fl(l+M 2 ) 2 * 

Example 2. By definition, a helix is a space curve whose tangent 
makes a constant angle with a fixed direction. Show that any single 
one of the following is a necessary and sufficient condition that a curve 
be a helix: (i) the principal normal is orthogonal to a fixed direction, 
(ii) k/t is a constant, (iii) r"xr"y"=0. 

(i) By the definition a helix is a curve such that 

T.e = cos a 


where e is a fixed unit vector and a is a constant angle. Differentiation 
yields T / .e = 0. Application of the first of the Frenet formulae 
followed by division by k yields N.e = 0. That is, the principal normal 
is orthogonal to the fixed direction e. Conversely, given N. e = 0, we 
have 

Jr (T.e) = T\e = *N.e = 0 


and so T.e is constant as required. 
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(ii) Differentiation of N. e = 0 yields ( — kT +rB). e = 0. Since N is 
orthogonal to e, both T and B are coplanar with e. Hence T. e = cos a 
and B.e = sin a if the relative positions are as depicted in Fig. 23. (If 



B is in the opposite direction, B. e = — sin a and the ensuing calcula¬ 
tions require slight modification.) Thus we have — k cos a + r sin a = 0 
and so /c/r=tan a = constant as required. 

Conversely, given k[t is constant, we can denote the constant by 
tan a to give k cos a—r sin a=0. Then we have 

Jr (T cos a+B sin a) = T' cos a + B' sin a 


= (k cos a—r sin «)N = 0, 

and so T cos a+B sin a = c, where c is a constant vector. Accordingly 
T. c = cos a and the curve is a helix. 

(iii) We have the relations 

r" = k N, 

T m = -^T + k'N + zctB, 

r"" = - 3 kk'T + ( k " - K 3 - kt2)N+(2k't+kt')B, 

and so 

^m ., 3 / / / \ 5 d l'r\ 

t .r XT’ = K % K T-KT) = K *‘-[-y 
Thus t/k = constant and r". r m x r ""= 0 each imply the other as required. 
Exercises 

L For the circular helix t=a cos u i+a sin w j+ bilk, show that 
k = a](a 2 +b 2 ) and r = b](a 2 +b 2 ). 

2. For the curve show that 

K = -r = V2 l(e u +e~ u ) 2 . 
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3. If a curve lie on a sphere show that (k 7 t#c2 ) = t / k * (Hint: The 

equation of the sphere with centre at c and radius a is (r—c). (r—c) = a 2 . 
Successive differentiation yields 


(r-c).T=0, (r—c).N=—1 /k, (r-c).B = k'/k 2 t 

and (r—c).N = 


T ds K Z T 


28 Surfaces and Normals 

Consider the vector function r(u, ©) of the two independent scalar 
variables u and ©.To each pair of values of u and © there corresponds 
a vector whose initial point is at the origin; the other end of the vector 
determines a point in space. As u and © vary this point traces out a 
surface. To avoid having the surface degenerate to a curve, the case 

when — x ^ = 0 is excluded. 
du dv 


All points on the surface for which © has the constant value © 0 lie 
on the curve r(«, © 0 ) with parameter u. This curve along which u 
varies is called a w-curve. Similarly all points on the surface for which 
u has the constant value u 0 lie on the ©-curve r(w 0 , ©). Collectively, 
the M-curves and ©-curves are called the parametric curves. 

The vector r u = drjdu is a tangent vector to the parametric curve 
© = constant whilst the vector r v — dr/d© is a tangent vector to the 
parametric curve u = constant. 

The equations u = u(t), © = ©(f) determine the curve r(«(f), ©(f)) 
with parameter t on the surface. It has the tangent vector 


dr 

dt 


dr du dr dv du 


dv 


du dt dv dt dt dt * v ‘ 


That is, at a point the tangent vectors to all curves on the surface 
passing through the point lie in the plane formed by r u and r v . This 
plane is called the tangent plane to the surface at the point. 

We define the normal at a point on a surface to be the straight line 
orthogonal to the tangent plane at that point. Let n denote the unit 
vector along the normal and be such that the vectors r u , r v and n form a 
right-handed system. The vectors r u and r v are both orthogonal to n 
since they lie in the tangent plane. It follows that 


r u xr v 
r u x r v \ 


n — 
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Example. Obtain the unit normal vector to the sphere 

r = a sin u cos v i + a sin u sin v j + a cos u k. 

(This vector equation does represent a sphere of radius a since the 
components x , y and z of r satisfy the relation x 2 +y 2 + z 2 — a 2 . 
The parameter u represents the colatitude (that is, the complement of the 
latitude) whilst v measures the longitude. The ranges of the parameters 
are 0 ^ u ^ tt and 0 ^ v < 2tt, The w-curves and ^-curves are respectively 
the meridians and the parallels of latitude.) 

We have 

r u = a cos u cos v i + a cos u sin v j — a sin u k, 
r v = —a sin u sin v i + a sin u cos v j, 

and so 

r u xr„ = a 2 sin 2 u cos vi + a 2 sin 2 u sin vj + a 2 sin u cos u k, 

|r u x r^l = a 2 sin m, 

which yields 

n = sin u cos v i + sin u sin vj + cos u k = r /«. 

Thus the positive direction of the unit normal is outwards from the 
sphere. 

Exercises 

1. Show that r=u cos v x + u sin v}+f(u)k represents a surface of 
revolution. What are the parametric curves ? Obtain the unit normal vector 
to the surface. 

2. Show that r =u cos v i -\-u sin vj + f(v) k represents the conoid, which 
is the surface obtained by revolving a straight line about an axis per¬ 
pendicular to it whilst simultaneously moving the straight line along the 
axis. What are the parametric curves ? Obtain the unit normal vector to the 
surface. 

29 Length of Arc on a Surface 

The distance ds between the two points r and r+dr on the surface 
r(u , v) is given by 

ds 2 = doc 2 + dy 2 + dz 2 
= dr.dr 

= (r u du+r v dv). (r u du + r v dv) 

= E du 2 + 2F dudv-\-G dv 2 y 

where 

E = r u .r u , F = r u .r v , G = r v .r v . 




4] LENGTH OF ARC ON A SURFACE 45 

E du 2 + 2 Fdudv + G dv 2 is called the first fundamental form of 
the surface r(w, v). 

In particular, the distance d u s along the M-curve for which dv = Q is 
given by 

d u s = \/E du. 

Similarly, the distance d v s along the w-curve is given by 

dyS = \/G dv. 

If every M-curve is orthogonal to every ®-curve, we say that the 
parametric curves form an orthogonal system. Hence the necessary 
condition that the parametric curves form an orthogonal system is 
that the tangent vectors r u and r v at each point be orthogonal; that 
is, if and only if F=r u .r v = 0. 

The scalar product —dr.dn plays an important role in the differen¬ 
tial geometry of surfaces and is called the second fundamental form. 

We have 


where 


— dr. dn = — (r u du + r v dv). (n u du + n v dv) 
= L du 2 + 2M du dv + N dv 2 , 


L — -r u .n u , 2M = - (r u . n„ + r v . n u ), iV = -r v .iv 

Since n is orthogonal to r u and r v we have r u . n = r v . n = 0. Differentia¬ 
tion yields 

^u • 4" r uu • n = 0, t u . n„ + r uv . n 0, 

r v . n u 4" r uv . n 0, t v . n^ 4" t vv . n 0. 

It follows that r u .n„ = r v .n u = —r uv .n and so we have 

■E = = \?v?v*uu\l\*u x 

M = T uv .n = [r u r v r uv ]l\r u xr v \ ) 

N = r vv .n = [r u r v r vv ]l\r u x r v \. 


Exercises 

\ 1. Obtain the first fundamental form of the sphere, surface of revolution 
and conoid. 

V2. Show that |r u xr„| — \/(EG—F 2 ). 

3. Obtain the second fundamental form of the sphere, surface of revolu¬ 
tion and conoid. 
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30 Scalar and Vector Element of Area 
Consider (Fig. 24) the four neighbouring points A, B , C, D on the 



surface such that the respective values of the parameters at these 
points are (u, v), (u + du, v), (u + du,v + dv) and (u, v + dv). The 
distances AB and AD measured along the parametric curves are 
respectively d u s=\/E du and d v s=\/G dv. The scalar element of 
area dS contained by ABCD is approximately a parallelogram and 
so dS=sin co dyS d^ where a> is the angle between the tangents to 
the parametric curves through A. Since r u and r v are tangent vectors 
to these parametric curves, we have 

|r u xr w | 
sin oi = \— 

| r tt| | r «| 

Consequently, we see that 

dS = |r u x r„[ du dv. 

The vector element of area dS is defined by 
dS = r„ x r„ du dv 


and 


so 


dS = n dS. 


(30.1) 

(30.2) 

(30.3) 






CHAPTER 5 
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31 Riemann Integral 

Let f(x) be a single-valued function of x defined in the interval 
a = x 0 ^x^x m =b. Divide this interval into m sub-intervals by the 
points * 0 , x lt —, x m where x 0 <x 1 <x 2 < — < x m - i< x m- Let us 
write Ax r =x r —x r _ 1 , choose £ r such that A: r _ 1 < € r ^x r and form the 
sum 

m 

I* = tttQAXr- 

r = l 

Let m tend to infinity in such a way that each Ax r tends to zero. If the 
limit of I m exists and is independent of the mode of subdivision, this 
limit is called the Riemann integral of f(x) from a to b and is 
denoted by 

f(x)dx. 

Ja 

It can be proved that this integral exists when f(x) is continuous 
in the interval {a, b). However, the integral may exist in cases when 
the function is not continuous. 


32 Line Integral 

The vector element of arc length ds is defined by 

ds = T ds, 

where T is the unit tangent vector drjds along the curve r($). Let T 
denote the portion AB of the regular curve r($) corresponding to the 
range s A of the arc distance s. Consider the vectors a(f) defined 

at all points of T. Then the line integral J* a.ds taken over T is 

defined to be the Riemann integral 


I a .ds — j B a.T ds. 
Jr Jsji 
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That is, the line integral of a vector is the Riemann integral along the 
curve of the tangential component of the vector. 

Let 

a = a^x, y, *)i + a^x, y, *)j + a k (x, y, *)k 

and 

r = #($)i + y($)j + #($)k 

referred to the fundamental system i, j and k. Then the line integral is 



where u A and u B are the parameter values of u corresponding to A 
and B. 

If r is a closed curve it is customary to denote the integral by 


o 

J J 


a 


ds. 


J. 

f^j/+sj + 


Example. Evaluate the line integral of^/+#j + ak over the curves 

(i) the circle x 2 +y 2 —1 = #= 0 , (ii) the triangle in the xy-plane with 
vertices at ( 0 , 0 ), ( 2 , 0 ) and ( 2 , 1 ), (iii) the skew quadrilateral with 
vertices at ( 0 , 0 , 0 ), ( 1 , 0 , 0 ), ( 1 , 1 , 0 ) and ( 1 , 1 , 1 ). 

The line integral is 


f / dx 


dy dz\ 
+ Z Tu +X d~u, 


du. 


(i) Over the circle traversed in the counterclockwise direction we 
have x =cos u, y = sin u, z =0 where 0 ^ u ^ 2tt and so 

r2n 

I z= — J sin 2 u du = —7 t. 

C dx 

(ii) Since »=0, we have /= y y du. From (0, 0) to (2,0) we 

have y = 0 and the contribution to the line integral is zero. From (2, 0) 
to (2, 1) we have x =2 and the contribution is again zero. From (2, 1) to 
( 0 , 0 ) we have x=—2u, y——u for the range — l^w^O and the 


contribution 


. . f° 

ion is 


2 udu— — 1. Thus 1= — 1. 


(iii) From (0, 0, 0) to (1, 0, 0) we have y = z= 0. From (1, 0, 0) to 
( 1 , 1 , 0 ) we have x=l t #=0 and so the contributions to the line 
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integral are zero. From (1, 1, 0) to (1, 1, 1) we have x=y= 1, z—u for 
the range 1 and the contribution to the line integral is 



1 . 


From (1,1,1) to (0,0,0) we have x=y = z=—u for the range 

— 1 <&<0 and the contribution is 3 J* udu— — 3/2. Hence the 

value of I over the skew quadrilateral is — 

This example shows that the value of a line integral usually depends 
on the choice of path. 


Exercises 

1. Calculate the line integral of yzi+zsq+xyk over the curves (i) the 
straight line joining the points (6, 0, 0) and (6, 0, 2r rc), (ii) the circular 
helix b cos ui + b sin u j+cwk which joins the above points. 

2. Calculate the line integral of (y 2 + z 2 )i +(# 2 4- x 2 )j + (x 2 +y 2 )k from 
(0, 0, 0) to (1,1,1) over the curves (i) the straight line joining the points, 
(ii) the three straight lines which link these points via the points (1, 0, 0) 
and (1, 1, 0), (iii) the curve wi+M 2 j4-M 3 k. 


33 Vector Line Integral 

With the notation of the previous section, divide T into m sub¬ 
intervals by the points with parameters s A =s 0 > s lf .. s m = s B . Denote 
the arc length of the sub-interval ($ r _i, $ r ) by As r — (s r — ^ r _ A ) and 
choose a r such that s r ^ x ^ oy < s r . Then form the vector sum 


I m — 2 a (°v) As r . 


r = 1 


Let m tend to infinity in such a way that each As r tends to zero. If the 
limit of I m exists and is independent of the mode of subdivision, this 
limit is called the vector integral of a over the curve T and is 
denoted by 


J* a ds. 

Further, the vector line integral J a x ds is defined by 
J axds = J axT ds. 


* 
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Example. Calculate (j) a x ds, where a=cos u i + sin u} + e u kand T 
is the circle p cos u i +p sin uj. We have s=pu and so 
r = p cos (s/p) i + p sin (s/p) j. 

Hence 

T = —sin (s/p) i + cos (s/p) j = —sin u i+cos w j 

and so 

a x T = — cos ue u i — sin ue u j+k. 

Accordingly the integral is 

J *2 n 

( — cos ue u i—sin ue u j + k) du = {(1 — e 2jl )(i +j) + 47 rk}. 

o 

Exercises 

1. Calculate (j) a ds for the vector a and the curve r of the above 
example. 

2. Show that <j) ax ds, where a=xi and r is the portion of the circular 
helix b cos u i+b sin u j + cwk between (—b, 0, ttc) and (b, 0, ttc), is 

%n 2 c 2 \ + 2bck. 


34 Double Integral 

Let f(x,y) be a single-valued function of the two independent 
variables x , y defined over a finite region S of the xy-plane. Divide 
the region (Fig. 25) into m sub-regions of areas AS X , AS 2 , .. AS m , 
(this is usually achieved by two families of curves), and choose any 



Fig. 25 
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point (£„ r] r ) inside or on the boundary of the sub-region AS r . Then 
form the sum 


4 = tm,Vr)JS r . 


r = l 


Let m tend to infinity in such a way that each AS r shrinks to a point. 
If the limit of I m exists and is independent of the mode of sub¬ 
division, this limit is called the double integral of f(x, y ) over the 
region S and is denoted by 

£/(*» y) dS. 


It can be proved that the double integral exists if f{x, y) is con¬ 
tinuous over the region S. 

The value of a double integral is computed by a suitable choice of 
sub-regions. For example, if the AS r correspond to the rectangles 
formed by straight lines parallel to the coordinate axes (Fig. 26) we 
write dS=dx dy. 



It can then be shown that 

£/(*> y)dS = j |J ^ f(x, y) dy J dx 

C d f r x I 
= i f{x,y)dx\dy i 
JC Uli(!/) J 

where y=y x {x) and y=y 2 (x) are the equations of the arcs ACB and 
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ADB whilst x=x x (y) and x=x 2 (y) are the equations of the arcs 
CAD and CBD. Further, a, b are the abscissae of the points A, B 
where the curve has vertical tangents whilst c , d are the ordinates of 
the points C, D where the curve has horizontal tangents. 

The regions considered in Figs. 25 and 26 are such that the 
boundaries are cut in at most two points by parallels to the axes. The 
definition of double integral can be extended to more complicated 
regions if they can be dissected into a finite number of sub-regions 
with the property just mentioned. 

Now let the vector a(x, y) be defined at all points of the region S. 
Then we can form the vector sum 

m 

*m = 2 ASr ' 

£r=l 

If the limit of I m exists as m tends to infinity and is independent of the 
mode of subdivision, this limit is called the vector double integral 
of a(x y y) over the region S and is denoted by 

J a(x, y) dS . 

35 Surface Integral 

Let a(w, ©) be a single-valued vector function of u and v defined 
over a region S of the surface r(«, v). Divide the region (Fig. 27) 
into m sub-regions of areas AS ly AS 2 ,..., AS m and choose any point 


| n (jUr.t'r) 

I 

I 



Fig. 27 
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(ti r , v r ) inside or on the boundary of the sub-region AS r . Then form 
the sum 


Ira v r) • v r) ^ 

r = l 

where n is the unit normal vector. Let m tend to infinity in such a 
way that each AS r shrinks to a point. If the limit of I m exists and is 
independent of the mode of subdivision, it is called the surface 
integral of a (u f v) over the region S of the surface r(u, v) and is 
denoted by 

I a.n dS = (* a.rfS, 

Js Js 


where dS = ndS. 

The surface integral can be evaluated as a double integral by 
interpreting u, v as rectangular Cartesian coordinates on a suitable 
region S of the w^-plane. In addition (see equation (30.1)) we replace 
dS by |r u x r„| du dv. Accordingly, 


£ a .dS = J_ a.n|r u x r„| dS t (dS = du dv). 

In particular, the area of any region S of the surface is given by 

L ds= L n - ds= S-s lw[dS - 

Tn | m 

2 ar I ^ 2 l^l an< ^ |a.n| ^ J a | it follows from the surface 


From 

Ir^l | r=l 

integral definition that 


|a.n| dS** J |a| dS. (35.1) 

A surface which encloses a finite volume will be called a closed 
surface. For example, a sphere is a closed surface but a hemisphere 
without its diametral plane is not a closed surface. For closed surfaces 
it is conventional to arrange that the positive direction of the normal 
is outwards from the surface. If S is a closed surface it is customary to 

denote the surface integral by (j) a . dS. 
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Three types of vector surface integrals are defined in terms of 
the corresponding vector double integrals by 


f. a 

I ax 


JS 



and 


fdS = \ fn |r u x r„| dS , 

Js Js 

where f(u, v) is a scalar function defined over S. 

Example. Evaluate J, a dS, J a.JS and J axJS for the vector 

a=cos u cos v i + cos u sin v j — sin u k over the octant of the sphere 
r= p sin u cos v i + p sin u sin vj + p cos u k corresponding to 0 < m < 7r/2 
and 0 ^7r/2. 

We can readily calculate that 


|r u x r„j = p 2 sin w, 

n = sin u cos vi + sin u sin vj + cos u k, 
a.n = 0 and axn = sin v i— cos vy 

Hence we have J a. JS = 0 and 


J a dS — p 2 J_ (cos u cos vi + cos u sin vj — sin u k) sin u dS , 

J a x JS = p 2 j_ (sin v i — cos v j) sin u dS, 

where S is the rectangle corresponding to 0 < m < 7t/2 and 0 < sy < tt/2 
in the wsy-plane. 

Simple calculations with dS = du dv now yield 
| adS = ip 2 (i+j-^) 


and 


J ax JS = p 2 (i-j). 


5] 
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1. Evaluate a dS, a . dS and a x dS for the vector 

Js Js Js 

a = (j'+*)i+(*+^j + (*4^y)k, 

where S is the surface of the cube bounded by a:=0,y=0, z=0, x— l,y== 1 
and z= 1. 

2 . Evaluate J xyz dS over an octant of a unit sphere whose centre is at 
the origin. 

36 Volume Integral 

Let f(x, y, z) be a single-valued function of x, y and z defined 
throughout a three-dimensional region V. Divide V into m sub- 
regions of volumes A V lf A V 2f ..., A V m . Choose any point (| r , -q rt £ r ) 
inside or on the boundary of the region A V r . Then form the sum 

r-1 

If the limit of I m exists as m tends to infinity in such a way that each 
A V r shrinks to a point and if the limit is independent of the mode of 
subdivision, the limit is called the volume integral of f(x, y, z) 
over the region V and is denoted by 

jf(x,y, z) dV. 

If the subdivision is made by planes parallel to the coordinate 
planes, the volume integral can be evaluated as an iterated integral 
by a generalization of the method outlined in Section 34. 

Corresponding to the vector function a(x, y, z), defined in the 
region F, we can form the vector sum 

r=1 

and so in the usual way by a limiting process obtain the vector 
volume integral 

Exercise 

1. Evaluate j a dV for the vector xi+yj+zk, where V is the region 
bounded by the surfaces #=0, y=0 , y = 6, z=4 and z=x?. 

5—V.A. 


CHAPTER 6 


Gradient of a Scalar 

Function 


37 Directional Derivative 

In this and the following chapters the function /(r) = f(x, y, z) will 
denote a single-valued function of the rectangular Cartesian co¬ 
ordinates x, y and z defined at all points of a region of space. We some¬ 
times say that a scalar field f has been defined in a region. 

We define the directional derivative D a f( r) of/( r) in the direction 
given by the vector a by 


DJ(r) = lim 

h-*0 


/(r + ha) — /(r) 
h 


provided that this limit exists. 

Let P and Q be the neighbouring points (#, y, z) and 

OJ* 

(x+Ax, y+Ay, z+Az) at a distance As apart. We use the notation — 

for the directional derivative of /(r) along the direction PQ. In this 
case h=As and a=Axi+Ay$+Az)s. and so 

^ = iim = lim f( x + Ax >y+ A y’ *) . (37.!) 

8s As -*0 As As -*0 

We have 


A f = (l +ei ) Ax+ % +e *) Ay+ (s +£3 ) Az> 

where e u e 2 and e 3 tend to zero as Ax, Ay and Az tend to zero. Divide 
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by As and proceed to the limit as As tends to zero with the result that 


df _ df dx df dy df dz 
8 s dx ds^dy ds^dz ds 


(37.2) 


It is clear from the definition of directional derivative that the 

partial derivatives and ^ are the directional derivatives in the 

dx dy dz 

directions of the coordinate axes. 

Simple calculations yield 

!(«/+« = a |+A g, 


where a and b are constants. 

The differential df is defined by 

dx dy 

and so 




K 

dz 


df = Qds. 

UJ ds c «. 


Exercise 

1. If/ =«: 3 +y 3 +a: 3 , find the directional derivative of / at the point 
(1, —1, 2) in the direction of (i) the vector i+2j+k, (ii) the vector j — k, 
(iii) in the direction of the normal to the surface ^c 3 +y 3 +« 3 =B. 

2. Find the directional derivative of |A| 2 in the direction a. 


38 Gradient of a Scalar Function 

The equation/( jc, y, z) = k represents a surface. If we vary the value 
of k we obtain a family of surfaces, one of which passes through any 
particular point. Let n be the unit normal vector to the surface. The 
directional derivative of f(x, y, z) in the direction of the normal to the 
surface is called the normal derivative and denoted by dfjdn. We 
define the gradient of f(x, y, z), denoted by grad /, to be the vector 

grad/= l n ' (38-1) 

Note carefully that this definition is invariant in the sense that grad / 
is independent of the choice of the basis. 
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It follows from the definition that grad c=0 if c is a constant. 

Let PN of length An be drawn (Fig. 28) along the normal to 



f(x, y, z) = constant at P and PQ of length As be drawn along a 
direction making an angle <p with PN such that N and Q are both on 
the same neighbouring surface of the family. Accordingly, NQ is 
approximately orthogonal to NP and so we have AnxAs cos cp. 

From the definition (37.1) we have that the directional derivative 
8f/8s along PQ is given by 


d l 

ds 


8f 

£ 008 ’’' 


Hence the maximum value of the directional derivative is attained in 
the direction of the normal. 

Let T denote the unit vector in the direction PQ. Then 


and so 


cos <p = T.n 


d l 

8s 



d l 

8n 


n.T = grad/. T. 


That is, the directional derivative in any direction is the component 
in that direction of the gradient. 

Further, we have 

df =if $ ds = grad/.T ds = grad/.ds, 
where ds = T ds. 


Let —> — and be the directional derivatives in the directions 

8s i 8 s 2 cs 3 


of the three mutually orthogonal unit vectors eu e 2 and e 3 . Then 

V 

8si 


grad/.e x = grad/.e 2 = JC grad/.e 3 = JC 
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since we have just seen that the directional derivative is the com¬ 
ponent of the gradient. Hence 


, , df df df 

grad/= e£ ei+ ^ ea+ ^ e3 - 


In particular, if we identify e 1} e 2 and e 3 with the fundamental 
system i, j and k, the corresponding directional derivatives are the 

partial derivatives ^ and Accordingly, we have 
dx dy dz 


, , df df . df , 

«" d /-S 1 + * J+ S k - 


From the relation gj=grad/.T and the properties of directional 
derivatives it follows that 


grad (af+ bg).T = ~ (af+ bg) 


,dg 

= a ds +b i 


where a and b are constants. 
Further, we have 


= a grad/. T+b grad g.T 
= (a grad f+ b grad g ). T, 


grad (fg).T = |(&) 


,dg df 
~ ^ds +g ds 


= / grad g . T +g grad/. T 
= (/ grad £+£ grad/).T. 

Since T is an arbitrary unit vector, we have 

grad (af+bg) = a grad f+b grad g, 
grad {fg) = / grad g +g grad/. 

Exercises 

1. Calculate the normal derivatives at (—1,1,1) of the functions 
(i) yz+zx+xy, (ii) xyz, (iii) x 3 —y 2 +z. 
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2. f Show that gradr=r/r and grad r n = nr n ~' 2 t. (Hint: Use a funda¬ 
mental system of basis vectors.) 

3. Show that grad (flg) = (g grad/-/grad g)/g 2 . 

4. If c is a constant vector, show that 

(i) gradc.r=c, 

(ii) grad |cxr| n =«|cxr| n_2 c x(r xc), 

(iii) c.grad ^c.grad={3(c.r) 2 — c 2 r 2 }/r 5 . 

5. Show that grad f(r) =/'(r)r/r. 

39 Irrotational Vector 

If a vector a(x, y, z) is defined at all points of a region, we say that 
a is a vector field in the region. 

The circulation of the vector field a round a closed circuit T is 

defined to be the value of the line integral (j) a.ds. The vector field a 

is said to be irrotational if its circulation round all circuits is zero. 

A region R is said to be connected if any two points of the region 
can be joined by an arc lying completely in the region. Figures 29(a) 
and ( b ) depict connected regions, but the region in Fig. 29(c) is not 
connected. 



t The vector r in this text will always denote the position vector relative to 
some origin. 
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A region is said to be simply-connected if every closed curve in 
the region can be shrunk continuously to a point in the region. The 
region R of Fig. 29(c) is simply-connected, whilst that of Fig. 29(6) 
is not simply-connected as the curve T cannot be shrunk continuously 
to a point without leaving the region. 

A necessary and sufficient condition that the continuous vector 
field a be irrotational in a connected region is that a single-valued 
function f(x, y, z) exists for which a = grad /. 

To prove the sufficiency, we note that 



since/is single-valued. 

To prove the necessity, let (j) a. ds — 0 and consider the case when 

the circuit T is the curve LMPNL in Fig. 30. The line integral round 
this closed circuit is zero and so we have 


J a.ds = I a.ds. 

LMP JLNP 



Choose L to be a fixed point. It follows that / a. ds taken along any 
path from L to P inside the region is a function of the position P 
and so we may write 

m = f 

Jlp 

The value of this function at a neighbouring point Q is given by 
f(Q) = f a.Js+ f a.ds. 

Jlp Jpq 

Hence the change Af in / from P to Q is 
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Select the straight line path from P to Q. By continuity, we have that 
along PQ 

a = a(P) + e 

where e tends to zero as Q tends to P. Accordingly 

Af=a(P).ds + f e.ds. 

Jpq 

The integral term is of higher order than ds and so in the limit we 
have 

df = a.ds. 

But df= grad f.ds. Thus 


(a— grad/).z/s = 0. 

This equation is true for arbitrary ds and so we have the required 
result that 

a = grad /. 

Example. Obtain the conditions that the vector 
x a y b z c (x l i +y m j + # n k) 

be irrotational. 

By the theorem just proved there must be a scalar function f(x, y, #) 
such that 


— /wCt + Z*. .6 f — Mtx.ib'i'fnnfC f ib«yC + ti 

gx -x- yz, gy-xy z, & *- . 

Integration yields the following three possibilities for /: 

x a+l + lyb z c 


a + l+l 


x ayb+m +lgc 

b V m+1 + *<*»*>» 


x a y b z c+n+ i 
c+n + 1 


+ 'l'3( x >y)> 


where ijj lt ifi 2 , */t 3 are suitably chosen functions of the variables indicated. 
The three values for / can be reconciled if 

(i) l=m = n— — 1, a = b = c. The corresponding scalar function is 
/= (xyz) a /a, 

(ii) a = b = c= 0 and the corresponding scalar function is 

xjl +1 v m + 1 ~n +1 

f =--1- 

l+l m +1 n+V 
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Exercises 

1. Show that (i) <|> t.ds=0, (ii) <|> c.ds — 0 for any constant vector c. 

2. It is given that P, Q and R are functions of x, y and z with continuous 
partial derivatives of the first order inside a simply-connected region which 
contains the path P joining (# 0 , y Q , z 0 ) to (*, y, z) by the straight lines 
connecting the points via (x, y 0 , z 0 ) and (x, y, z 0 ). Show that 

a = Pi+0j+Pk 

is irrotational if and only if 

8R 8Q 8P _8R 8Q 8P 

8y 8z * 8z 8x an 8x 8y 


(Hint: define A(x,y, z) = J a.ds and show that a=grad A.) 

3. Determine which, if any, of the following vector fields are irrotational: 
(i) (X.r)fx, (ii) (X.fi)r, (iii) (X.r)p.+(p..r)X, (iv) (X.r)p.-(p..r)X, where X 
and (X are constant vectors such that X x p.^0. In appropriate cases, obtain 
the scalar field whose gradient is the irrotational field of vectors. 


40 Integral Definition of Gradient 

Let S be a closed surface which contains a volume V and define 
the vector a from the scalar field f(x,y, z) by the equation 


a = 



provided that the limit exists. 

Choose S to be a cylinder (Fig. 31) with its generators parallel to 
the fixed unit vector c. Then 


a.c 


= lim C) /c 
v-*o V J S J 


.dS. 



Fig. 31 
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On the curved surface of the cylinder c is orthogonal to the unit 
normal vector n and so c. dS = c. n dS= 0. Thus the total contribution 
to the surface integral comes from the base and top of the cylinder. 

Let the cylinder be of height As and cross-sectional area AA. On 
the base c. dS = c. ( — c dS)— — dS whilst on the top 


c.d'S = c.c dS = dS. 


By the mean-value theorems, the contributions of 
base and top are 



.d'S to the 


~(f+ei)AA and [f+^As + e^ AA, 
where e x and e 2 tend to zero as AA tends to zero. Thus 


™0 As A A ~ di 

AA-*0 


But the component of grad / in the direction c is ^ and so 

a = grad/. 


That is, grad / can be defined as the limit of an integral. Since the 
unit vector c is arbitrary, it follows from this definition that grad/is 
independent of the choice of the basis vectors. 
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41 Divergence of a Vector 

Let a be a vector field defined inside a closed surface S which 
encloses a volume V. We define a scalar field associated with a, called 
the divergence of the vector a and written div a, by the relation 

(j) a.dS 

div a = lim - — (41.1) 

v-*o r 

provided that the limit exists. It is clear from the definition that div a 
is independent of the choice of the basis vectors. 

If A, /x and c are constants, we have from the definition that 


and 


div (Aa+ftb) = A div a+/z div b 
div(/c) = lim (j) /c.dsj- = c.lim (j) fdS 


That is, 


div (/c) = c.grad/. 


Substitute /— 1 and we obtain div c = 0. 

Let a = a t i + a } j + c fc k with respect to the fundamental system i, j, k. 
Then we have 

div a = div (af) + div («J) + div (« fc k) 

= i. grad a t + j. grad a t +k. grad a k . 

But i.grad a t is the directional derivative along the ac-axis, etc. Thus 
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The Laplacianf Af is defined by 

Af = div grad / 

and so in a Cartesian coordinate system we have 

A f= &LJl&L. 

J dx 2 8y 2 8z 2 

Exercises 

1. Show that div r n r=(«+3)r n . 

2. If c is a constant vector, show that (i) div (c x r) = 0, (ii) div (c. rc) = c 2 , 
(iii) div {c x (r x c)} = 2c 2 . 

3. If p=x 2 i+y 2 j+z 2 k, q=yzi+zxj + xyk and u=xyzj, calculate 
(i) div p, (ii) div q, (iii) div u, (iv) div (qxu), (v) grad (p.q), 
(vi) grad div (q x p), (vii) A [pqu]. 

4. Calculate (i) A(ljr), (ii) dr 2 , (iii) Af{r), (iv) d[div (r/r 2 )]. 

42 Gauss’s Theorem 

Consider a finite closed surface S enclosing a volume V. Divide V 
into m sub-regions of volumes A V lt A V 2 ,..., A V m and choose any 
point (£ ry 7) ry £ r ) inside or on the boundary AS r of the sub-region A V r . 
From the divergence definition we have 

AV r div a(£ ry v „ Q = c£ a.dS+e r AV„ 

J AS r 

where e r tends to zero as AV r tends to zero. Addition of all these 
equations for r = 1 , 2, .. m yields 

m m /• m 

]?AV r div a(£ r , Vn Q = J I a.dS+ J e r AV r . 

r = 1 r= l J AS r r = 1 

The boundary of each A V r consists of a number of pieces which are 
either part of the boundary of S or part of the boundaries of two 
adjacent sub-regions. The total contribution from two adjacent 
boundary surfaces is zero since the outward normals have opposite 
directions over the common boundary surface. 

m 

Further, we have ^ e r A V r < V max e r and so this term tends to 

r = l 

t No confusion will ensue from the use of Af in two different ways. It will 
be clear from the text when A is the Laplacian operator and when it refers to an 
increment. 

Many authors use the notation V s instead of A for the Laplacian operator. 
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zero as A V r tends to zero and m tends to infinity. Thus in the limit 
we have 


J* divadT = <j) a.dS. 


(42.1) 


This result, by which a volume integral is transformed into a 
surface integral, is known as Gauss’s theorem or the divergence 
theorem. 

Since div (/c) = c.grad/, where c is an arbitrary constant vector, 
the substitution a=/c in Gauss’s theorem yields 

c. J grad/dT = c. <j) fdS. 

But c is an arbitrary vector and so we have 

J* grad fdV = <j) fdS. 

Choose/= 1 and we deduce that for a closed surface 

dS = 0. 


(42.2) 


<) 
J i 


(42.3) 


Let us apply Gauss’s theorem to the vector field a defined in the 
volume V contained between the closed surface S and a sphere 2 of 
radius R, centre at the origin, and sufficiently large to enclose S 
completely. We have 

divadF = <p a.dS+ u) a.dS, 

JV J S J 2 

where the normal in the surface integrals is directed outwards from V. 
From equation (35.1) we have 


, H ds < W 


% 

() 
J 2 


| j) a.dS ^ ^ 

if |aj ^ ejR 2 at all points on 2* But (j) dS is the area of the sphere of 


radius R and so 


\i 


a.dS 


^47 re. 


Now let us impose the condition on a that lim R 2 a=0. Then c 

R-*co 

can be made as small as we please and we obtain 
J* divadF = £ a.dS. 
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That is, Gauss’s theorem extends to the infinite volume outside a 
closed surface provided that lim R 2 a = 0. 


R-* co 


Example. Evaluate <j) (a 2 x 2 + b 2 y 2 + c 2 z 2 ) ~ y * dS over the surface of 

the ellipsoid ax 2 + by 2 + cz 2 = 1 . 

The normal to the ellipsoid is along the direction of 

grad (ax 2 + by 2 + cz 2 ) — 2(axi + by] + c#k). 

Thus the unit outward normal vector n to the ellipsoid is 

n = (axi + by] 4- czk)j(a 2 x 2 + b 2 y 2 + c 2 z 2 ) 112 . 

Accordingly, the surface integral is <j) r-^S, where r=#i+jj+;sk 

since ax 2 + by 2 + cz 2 = 1 on the surface. 

By Gauss’s theorem, the surface integral is equal to the volume 

integral J div r dV — 3J dV. The volume enclosed by the ellipsoid 

is 47rj[3\^(abc)] and so the original surface integral has the value 

47r/V(«^)* 

Exercises 

1. Evaluate <j) a.dS, where a =yzi + + *yk and S is the surface of 

the cube formed by #=0, x— 1, y= 0, y= 1, z=0 and z= 1. 

2. For a closed surface S containing a volume V, show that 


(i) £ r.dS=3V, (ii) i dV=j 


r.dS 


3. Evaluate <p ( ax 2 + by 2 + cz 2 ) dS , where S is the sphere x 2 +y 2 +z 2 = 1. 

4. Show that the value of j> p(^ 2 +p+^ dS taken over ^ surface 

of the ellipsoid x 2 la 2 +y 2 jb 2 +z 2 lc 2 =l , where p is the perpendicular from 
the origin to the tangent plane at (x, y, z), is ATrabc{a 2 +b 2 + c 2 )^. 

43 Divergence of the Product of a Scalar and a Vector 

Corresponding to the vector field a(ac, y, z) and the scalar field 
f{x, y, z), the divergence definition yields 


div (/a) = lim /a-dsj 
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Let the values of a and / at some fixed point inside the small 
region V be a 0 and/ 0 . Then we may write a=a 0 +Aa and f=f 0 +Af 
on the surface S. Hence 


div (fa) = lim (f 0 +Af)a.dS 
v-*o v J S 

= Mm y j/o (j>^ a . dS + j>^ (/-/ 0 )(a 0 +Aa ). j 
= limi|/ 0 (j) a.</S + a 0 . fdS 


* r ^ 

-/oa 0 .() dS+o AfAa.dSV. 
• s « s J 


The last integral on the right-hand side will vanish in the limit as it 
is of higher order than the other terms within the braces. Further, 

(f> dS=0 by equation (42.3). Hence 


div (/a) = Urn y^fo a.</S+a 0 .j> /</s\. 

That is, 

div (fa) — f div a + a. grad/. (43.1) 

This result can be obtained more easily in the special case when the 
basis consists of three mutually orthogonal fixed unit vectors. Let the 
basis vectors be e 2 and e 3 and let a =+ a 2 e 2 + a 3 e 3 . 

From div (/c) = c. grad/ in the case of the constant vector c, we 
have 


o 

div (/a) = J div (fa r e r ) 

r = 1 

3 

= 2 e r*8 rad (f a r) 


r = 1 
3 


= /2 e - S rad «r + a. grad/. 

r — 1 

Substitute f—1 and obtain 

3 

div a = 2 e r .grad a r 

and so 


div (/a) = / div a + a . grad /. 


Exercises 

1. Obtain f(r) if div {/(r)r}=0. 

2. Show that div grad (fg)=f div grad g+g div grad f+2 grad/.grad g. 
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Let a(x t y, z) be a vector field defined inside the region V enclosed 
by the surface S. We define a vector associated with a, called the curl 
of the vector a and denoted by curl a, by the relation 


<p axdS 

curia = — lim ■ s -y - (44.1) 

v-*o V 

provided that the limit exists. 

It follows immediately from the definition that 

curl (Aa+/xb) = A curl a+/x curl b 

if A and /x are constants. If c is a constant vector we have 


curl(/c) = — lim | /cx Jsj- = -cxlim <j) /Jsj-. 
That is, 

curl ( fc) = grad /x c, (44.2) 


from which we deduce that curl c=0. 

Let a=a t i + a i } + aji with respect to the fundamental system i, j 
and k. Then 


curl a = curl {af) + curl (%j) 4- curl {a k k) 

= grad a t x i 4- grad % x j 4- grad a k x k 


\&r dy^ dz ) \cfc dy^ 8z ) ^ 

■ (8 a k • ■ & a k • ■ & a k V\ « 

+ [dx 1+ 8y i+ 8z k j Xk 
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- (t-SHS-SMSf-tH <«> 

We deduce that 

curl grad/= 0 and div curl a = 0 

for arbitrary / and a. 

Exercises 

1. If a=zi~xj +yk, b=yi+zj—xk and /= xyz, calculate (i) curl a, 
(ii) curl b, (iii) curl grad /, (iv) div curl a, (v) curl (grad /x(axb)}, 
(vi) div grad/, (vii) curl curl (axb). 

2. Show that curl (r n r)=0. 

3. If a—g grad/, where/and g are scalar fields, show that a. curl a=0. 

4. If X and fx are constant vectors, show that 

(i) curl (X x r) =2 X, (ii) curl (fx.rX}=Xxfx, (iii) curl (X x (r xfx)} = jx x X, 
(iv) curl (r x (jx x r)}=3r x jx. 

5. By a method analogous to the proof of Gauss’s theorem show that 

curl a dV— — (j) ax dS. Deduce that J grad /dV = f dS. 

6. Show that j) rxdS = 0. 

7. For any scalar field/, show that 

(i) ^/grad/x dS=0, (ii) j> f(r)r x dS = 0. 

45 Curl of the Product of a Scalar and a Vector 

Corresponding to the vector field a(x,y, z) and the scalar field 
f( x , y, #), the curl definition yields 

curl (/a) = — lim j> fax dsj-. 

With the notation of Section 43 and by a similar method we obtain 

curl (/a) = -lim i|/ 0 ^axdS + a 0 x j) fdsj. 


t This expression for curl a with respect to a fundamental system can be 
written symbolically in the determinantal form 

i j k 

e_ d_ e_ 

iSx 8y 8z\ 

\ a t o,j a^t 


6—V.A. 
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That is, 

curl (/a) = /curl a + grad/x a. 

(45.1) 

Exercises 



1. Show that curl {/(r)r}—0. 

2. Calculate curl b in the cases when 

(i) b=e x+v+z (yzi+zty+xyk), (ii) b=xyz(e x i+e 3/ j+e 3 k). 



46 Divergence of a Vector Product 

Corresponding to the vector fields a(x, y, z) and b(#, y, z), the 
divergence definition yields 

div(axb) = lim axb.JS 

V —0 IT J s 

With the notation of Section 43 and by a similar method we obtain 


div(axb) = lim axb 0 .</S + j) a 0 xb.</sj-. 

We have 

axb 0 .JS = — b 0 .axJS and a 0 xb.JS = a 0 .bx dS. 
Consequently 

div(axb) = lim b 0 .^ ax dS + a 0 .(j) bx^sj- 


and so 


div (a x b) = b. curl a — a. curl b. 


(46.1) 


Exercises 

1. If c is a constant vector, show that (i) div |cx grad -^=0, 

(ii) div {c x/(r)r} = 0, (iii) div {r x grad/(r)}=0. 

2. If X and jx are constant vectors, show that (i) div {rx(Xxr)}=2r xX, 
(ii) div {(X x r) x (p x r)}=2X. ([x x r)—2(x. (X x r). 


47 The Operator a.grad 

Let a(x , y, z) and b(#, y, z) be vector fields defined inside the 
region V enclosed by the surface S. We define the vector associated 
with both a and b, denoted by (a.grad)b, by the relation 

(£ b(a.JS) 

(a.grad)b = lim ——-» (47.1) 

v-»o V 
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where a is to be regarded as constant over the surface S and provided 
that the limit exists. 

It follows immediately from the definition that 


(a.grad)(Ap+/*q) = A(a. grad)p+/x(a. grad)q, 

where A and [i are constants. Further, if c is a constant vector, we 
have 

(a.grad)(/c) = lim (|) /c(a.dS) = ca.lim (|) fdS, 
v-*o V J s v-»o V J s 


since a is constant on the surface of S. That is, 

(a.grad)(/c) = c grad/.a, 

from which we deduce that (a. grad)c = 0. 

Let a = a t i + aj j + a k 1z and b = b t i + b t j + ft fe k with respect to the 
fundamental system i, j and k. Then 

(a.grad)b = (a.grad)(fe i i) + (a.grad)(^j) + (a.grad)(fe fc k) 

= i grad b t . a+j grad fe^.a + k grad b k . a 

J8b, 8b,. 8b, t \ . . tx ... 

= l l +—j + — k j. (a t l + a i f j + a k k) + two similar terms 

/ 8b, 8b, 8b,\ . ... 

= ^ ~dx + aj ~3y + aic ‘dz) 1 + tw ° sim “ ar terms 


8b 8b 


8b 


~ 01 8x + a ‘ 8y + ak 8z 


Exercises 

1. Show that ax curl a=\ grad a 2 — (a.grad)a. 

2. Show that (a.grad) ^b.grad^={3a.rb.r—r 2 a.b}/r 5 if a and b are 
constant vectors. 

48 Gradient of a Scalar Product 

Corresponding to the vector fields a(x, y, z) and b(#, y, z) the 
gradient definition yields 

grad (a.b) = lim a.bdS. 
v-o V J s 
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With the notation of Section 43 and by a similar method we obtain 

grad (a.b) = lim <^(j) a.b 0 JS+ <j) a 0 .b JSj-. 

On application of the vector triple product formulae 
a 0 x(bx dS) = a 0 .dSb — a 0 .b JS, 
b 0 x(axJS) = b 0 .dSa — a.b 0 dS, 

we have 

grad(a.b)= limb 0 x<p axJS —aoX(p bxJS 
v-o V ( J s J s 

+ | a(b 0 .<JS) + (J> b(a„.<iS)}. 

Consequently 

grad (a. b) = b x curl a + a x curl b + (b. grad)a+(a. grad)b. (48.1) 

49 Curl of a Vector Product 

Corresponding to the vector fields a(x, y, z) and b(x, y, z), the 
definition of curl yields 

curl (a x b) = — lim (p (a x b) x dS. 
v-o V Js 

With the notation of Section 43 and by a similar method we obtain 

curl(axb)= — lim (a 0 x b) x JS + (j) (axb 0 )xJsj-. 

Expanding the vector triple products we have 

curl(axb) = —lim 77 /(p ba 0 .dS — <p a 0 b .dS 
v-o V U s J s 

+ <j) b 0 a.JS — <j) ab 0 .Jsj-, 

and so 

curl (a x b) = a div b - b div a + (b. grad)a - (a. grad)b. (49.1) 

Exercises 

1. If c is a constant vector, show that grad (c x r) = c. 

2. If c is a constant vector, show that 

curl x grad ^ + grad . grad ^ = 0. 
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50 Alternative Definition of Curl 

In the curl definition (44.1) choose V to be (Fig. 32) a right circular 
cylinder of cross-sectional area e and height /, where / 2 is small 
compared with e. Let the unit vector in the direction of the generators 
of the cylinder be c. Then 



Fig. 32 


The two flat ends of the cylinder make no contribution to the 
integral since c and n are parallel there. At a point P on the curved 
surface, c and n are mutually orthogonal unit vectors and so c x n=T, 
where T is the unit tangent vector to the plane section T P through P 
parallel to an end. By convention, n is the outward normal and so the 
direction of T P is determined by the fact that c, n and T form a 
right-handed system. 

We now evaluate the surface integral as a double integral. Let A 
be the length along a generator measured from the base of the cylinder 
to P , and s be the arc length along r, the base curve of the cylinder. 
We note that T P is obtained by translating T through a distance A 
parallel to c. 
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Since a depends on both s and A, whilst n is independent of A, we 
have 

<f) a.c x n dS = <|) a.cxndXds 
Js JrJ a=o 

= f r cxn [r a dX j ds. 

As / 2 /c tends to zero with e, Taylor’s expansion (Section 22) 
allows us to write 

a(s, A) = a(s, 0 )+iq(s, A) 

along a generator, where tq tends to zero with A. Accordingly, 


J a dX = /(a + jx), 


where p, is a vector function of s and / and tends to zero with /. 
Noting that V =e/, we finally obtain 


lim 

v-»o 

1 

el\ 

j) cxn./(a + p,) 

lim 

H 

p cxn.a ds 

£-»0 

«J 

* r 

lim 

1 j 
- < 

b a.ds. 

e—»0 

« J 

' r 


(50.1) 


Accordingly, the component at a point in a direction c of curl a is 
the limit of the ratio of the circulation of a around a plane curve r 
normal to c and containing the point to the area enclosed by r. 

The direction of ds in the circulation integral is such that c, the 
outward normal to r and ds form a right-handed system. 


Exercises 

1. Determine curl a in rectangular Cartesian coordinates by choosing T 
to be a rectangle in a coordinate plane. 

2. If a is an irrotational vector, show that curl a=0. Hence deduce that 
curl grad /= 0 for arbitrary /. 


51 Stokes's Theorem 

Consider a space curve r bounding an open surface S. Divide S 
into m sub-regions so small that they may be considered to be planar 
with areas AS Xi AS 2 ,..., AS m and choose any point (£ r , ^ r , £ r ) inside 
or on the boundary T r of AS r (Fig. 33). 
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Choose a positive sense of description for T. Then an orientation for 
each r r is determined by the conditions: 

(i) if r r and jT have an edge in common, this edge is described in 
the same direction along both boundaries, 

(ii) if r r and jT s (r ^ 5 ) have an edge in common, this edge is described 
in opposite directions. 

Let the unit normal at (£ r , r) r , £ r ) be n r with positive direction such 
that a right-handed corkscrew rotating in the direction of jT r will 
move in this direction. 



From equation (50.1) we have 

n r .curl a(£ r , rj r , t, r )AS r — (j) a. ds + e r AS n 

where e r tends to zero as AS r tends to zero. Addition of these equations 
for r = 1 , 2 ,. .., m yields 

m m n m 

y n r . curl a(£, v „ Q AS r = Y <b a .ds + V e r AS r . 

r = 1 r = 1" r = l 

The boundary jT r of each AS r consists of a number of pieces which 
are either part of the boundary T or part of the boundaries of two 
adjacent sub-regions. The contribution to the circulation from the 
two adjacent boundary curves is zero as they are traversed in opposite 
directions. 

m 

Further, we have ^ e r AS r < S max c r , where S is the total area of 

r = 1 

the surface, and so this term tends to zero as m tends to infinity in such 
a way that each AS r shrinks to a point. Hence in the limit wejiave 
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where n denotes the vector field of positive unit normals to the surface 
S. That is, 


J* curl a. dS = <j) a. ds. 


(51.1) 


This result by which a surface integral is transformed into a line 
integral is known as Stokes’s theorem. 


Exercises 

1. Verify Stokes’s theorem by evaluating both integrals for the vector 
yi+zj + xk and the surface S given by x 2 +y 2 +z 2 =l and z^O. 

2. Substitute a=/c, where c is a constant vector, in Stokes’s theorem 

and deduce that I grad/xdS= — j) fds. 

3. If S is a closed region of the ary-plane bounded by a simple closed 
curve C, deduce from Stokes’s theorem that 

!M- a 4) ds= L {Adx+Bdy) - 

This result is known as Green’s theorem in the plane. 


52 Surface Integral of the Curl of a Vector 

Consider a closed surface S (Fig. 34) and let the closed circuit T 
divide it into two regions S ± and S 2 . 



Apply Stokes’s theorem to both S ± and S 2 and the results are 


J curl a. dS — (j) a. ds 

I curl a. dS = (b a. ds. 
Js a J r 


and 
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The orientations of the normals to S x and S 2 differ, one surface 
having the normal in the outwards direction and the other in the 
inwards direction. Hence by subtraction we have 

<j> curl a.JS = 0 

for a closed surface S. 


53 Curl of the Gradient of a Vector 

Substitute a = grad /, where f(x, y, z) is a scalar field, in Stokes’s 
theorem to obtain 

j curl grad/. dS = <j) grad/.Js = 0 

in virtue of the fact that grad / is always an irrotational vector as 
shown in Section 39. However, S is any arbitrary surface and need 
not be closed. It follows that 

curl grad / = 0 

for any function /. 

If the vector field a satisfies the relation curl a=0, it follows from 
Stokes’s theorem that a is irrotational. In such a case we have shown 
(Section 39) that a scalar function f{x, y , z) exists such that a =grad/. 
It follows that curl a=0 implies the existence of a scalar function/ 
for which a=grad /. 

In this proof an appeal has been made to Stokes’s theorem. It 
follows that the region in which a is defined must be simply-connected. 

We now show by a counterexample that the vector field a may not 
be irrotational if curl a = 0 in a region which is not simply-connected. 
Consider the vector 


y x 

® O . O 14" O Q J* 

x 2 +y 2 x 2 +y 2J 

This vector a is not defined at the origin but is defined in the multiply- 
connected region between any two spheres with centres at the origin. 
In this region we readily verify that curl a = 0. 

The line-integral of a round any closed circuit r is given by 


- i, 
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Select r to be the circle x 2 +y 2 = p 2 , z = 0. The substitution 
x — p cos xfj, y = p sin ip yields 

J '2n 

dip = 27r ^ 0 

0 

and so a is not irrotational. 

Exercise 

1. If the vector field a is irrotational in the region V enclosed by the 
surface S and if a is along the surface normal on S, show that 


for any scalar field /. 


j ax grad fdV = 0 


54 Divergence of the Curl of a Vector 

Replace the vector a by curl a in Gauss’s theorem, equation (42.1), 
and the result is 


div curl a dV 


curl a. dS = 0 


since S is a closed surface. The volume integral is zero over any 
arbitrary region and so we have 

div curl a = 0 

for any vector a. 


55 Solenoidal Vectors 

A vector field a is said to be solenoidal if div a = 0. Since 

div curl a = 0 

for all vectors a, it is clear that curl a is always a solenoidal vector. 
We now establish by a method due to L. Brandf the converse result 
that corresponding to a solenoidal vector field a there exists a vector 
field b such that a — curl b. 



Fig. 35 

f Amer. Math. Monthly , 57, 161 (1950). 
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Consider the closed surface S formed by a cone with vertex at O 
and plane base (Fig. 35) determined by the closed circuit r. Since 
div a=0, Gauss’s theorem yields 


x 


a.dS = 0. 


Let us denote the curved surface by S lt the flat base by S 2 and the 
equation of r by r= r(s), where s measures arc distance along r and 
r has its initial point at O. The equation of S l is R=Ar(s) for the 
range A< 1. The vector element of area dS in the direction of the 
external normal to S x is given by 

gR dR , _ „ 

dS = -rr x —r- ds dX = AT x r as aA, 
cs OA 

where T is the unit tangent vector to r. Thus 


* 

a.dS = IT Aa(Ar).Txr<fcdA 

Si JJs i 


= <j> J Aa(Ar) <£\j .Txr ds 


= <j) rx Aa(Ar) dA| .ds. 

By application of Stokes’s theorem on the assumption that a=curl b 
we have 


1 a. dS = J curl b. dS = <p b. ds. 

s z Js z J r 

Accordingly 

/% /» 

where 

0 = y a.dS = y p.ds, 


p = b + r x J Aa(Ar) dA. 

Thus the vector field p is irrotational and a scalar function f(x, y, z) 


exists such that p = grad/. Hence 

b = grad / — r x J Aa(Ar) dA. 


(55.1) 


That is, we have obtained a vector b for which a=curl b. Note 
that b is indeterminate to the extent that we may add the gradient of 
any scalar function. We call b a vector potential of the solenoidal 
vector a. 
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Example. Show that the vector a = si+xj +yk is solenoidal and 
obtain a vector potential for it. 

Clearly diva=0 and so a is solenoidal. Application of equation 
(55.1) shows that a vector potential is 

— rxj A(A#i+Axj + Ayk) dX 

— ~~ 3(^1 +yj + #k) x (^i+xj +yk) 

= -^{(;sx —y 2 )i +{xy — z 2 )j + (yz— x 2 )k). 


Exercises 

1. Show that the vector grad/x grad g corresponding to the scalar 
fields / and g is solenoidal. 

2. Show that the following vectors are solenoidal and obtain vector 
potentials: (i) (y—z)i+(z— x)j + (x— y)k, (ii) r 3 cxr, where c is a constant 
vector, (iii) yzi + zoq + xyk, (iv) e x j. 

3. If the vector field a is both irrotational and solenoidal, show that 
grad (a. c) = curl (a x c) for any constant vector c. 

4. If the vector field a is solenoidal in the region V enclosed by the 
surface S and if the component of a along the surface normal vanishes, show 

that J a.grad fdV= 0 for any scalar field/. 

5. If the vector field a is solenoidal, show that 2J r.a dV=(£ r 2 a. dS. 
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56 Green’s Theorems 

In the relation div (/a)=/div a+a. grad/ let us write a=grad g, 
where g(x, y, z) is a scalar field. Then 

div (/grad g) = /div (grad g) + grad g. grad/ 

= fAg+grad f. grad g. 

Integrate over the region V contained by the closed surface S and we 
obtain 


f div (/ grad g) dV = f fAgdV+ f grad /. grad g dV. 

Jv Jv Jv 

The integral on the left-hand side of this equation is l fgradg.dS 
by Gauss’s theorem. However, 

/grad^.JS = / grad^.n dS = f^dS 
and so we obtain 


£/ff ^ = \y fAg dV+ jy 8 rad /-g rad ^ dV ' 

(56.1) 

We deduce that 


£/§£ dS = jjAfdV+j r (grad/) 2 dV 

(56.2) 

and 



(56.3) 


These relations are generally known as Green’s theorems. 
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57 Harmonic Functions 

A scalar function <p is said to be harmonic if it satisfies Laplace’s 
equation A<p — 0. 

Let <p x and <p 2 be two harmonic functions in the region V enclosed 
by the surface S such that 951 = 952 on Then f=<p x — <pz is harmonic 
in V and vanishes on S. Application of the equation (56.2) yields 

J (grad/) 2 dV = 0. 

Hence grad/=0 in V. That is, / is constant in V. But / is zero on S 
and so /= 0 in V. That is, 9 >i= 9 > 2 in V. 

Consequently, a harmonic function in a region bounded by a 
closed surface is uniquely determined by its values on the boundary. 

Substitute /= 1 and g = <p in equation (56.3) when <p is harmonic. 
The result is 

§M ds =°- {57A) 

That is, the surface integral of the normal derivative of a harmonic 
function is zero over the bounding surface of the region in which it is 
harmonic. 

Exercises 

1. If 95 is harmonic in the region V enclosed by the surface S and 9 ? is 
constant on S, show that <p is constant in V. 

2. If <p x and <p 2 are both harmonic in the region V enclosed by the 
surface S and their normal derivatives are equal on S, show that tp x —tp 2 is 
constant in V. 


58 Uniqueness Theorem 

We now prove that there can be at most one vector field in a region 
V contained by a closed surface S, if the divergence and curl are 
given everywhere in V and the component of the vector in the direction 
of the normal to the surface is given on S. 

Suppose a and b are two vectors satisfying the given conditions. 
Form the vector p = a — b and we have 


and 


div p = 0 and curl p = 0 in V 
0 on S. 


p.n = 
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Since curlp = 0, the vector p is irrotational and thus a scalar field 
f(x, y, z) exists such that 


Accordingly 


p = grad/. 


and 


Af = div grad/ = div p = 0 in V 


Application 


d l 

dn 


£ = grad/.n = p.n = 0 


of equation (56.2) yields 


on S. 


f (grad/) 2 dV = 0. 

The integrand is never negative and so we have (grad f) 2 — 0. That is, 
p = grad/=0. Hence a=b as required. 

Note carefully that we have not shown that a vector field exists 
satisfying the given conditions. We have merely shown that if one 
such vector field exists, then it is unique. 


59 Solid Angle 

Form the cone joining the point O to the points of a closed curve 
r. The solid angle of the cone is defined to be the surface area of the 
unit sphere, with centre at O, intercepted by the cone. 

Consider the vector 


11 1 

a = grad- = -^gradr = -^r, (r ^ 0), 


where r refers to O as initial point. Further, in virtue of equation 
(43.1) we have div (/b)=/divb+b. grad /and so for r#0 we have 


diva = 



—p div r—r.grad 




Apply Gauss’s theorem, equation (42.1), to the region (Fig. 36) of 
the cone bounded by the curved surface G formed by the generators, 
an open surface S whose boundary is r and the portion a of the small 
sphere of radius p intercepted by the cone. (We emphasize that p will 
not tend to zero. We choose p sufficiently small so that the surfaces S 
and a do not intersect.) The result is 



~.dS + 
r 



~.dS + 
r 



-.dS = 0. 
r 
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On G we have grad-= — i r and r is orthogonal to the surface 
T T° 

normal n. Accordingly, f grad-.dS=0. On a we have r=p and 

Jg r 

n= —rjpy the sign being negative because the normal is outwards 
from the surface, and so 



But the solid angle Q of the cone is given by J'jS' taken over the 
portion of the unit sphere intercepted by the cone. Hence J* dS=p 2 Q 
and so J grad ^-.dS = jQ. Thus we obtain 

Q = — J* grad^-.dS. 

As a corollary, it follows for a closed surface S that 


i 


grad — .dS = — Att if O is inside S, 
s r 


= 0 


if O is outside S. 


Exercise 

1. Show that the solid angle subtended at the origin by a closed curve r 
z(y dx—x dy) 

is given ~ 


by ±<|> 


i r {x? +y 2 )(x 2 +y 2 + z 2 ) 112 

60 Green’s Identity 

We apply equation (56.3) with S = ~ f where r is the distance from a 
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fixed point O, to the region V contained by the closed surface S. In 
virtue of the result Ag = div grad ^ = 0 established in the preceding 
section we have 


Consider the case when O lies inside the surface S. We exclude 
O by surrounding it with a small sphere a of radius e. Then we have 

where V is the region contained between a and S. 

On the surface a we have r=e and the external normal is directed 


towards O. Thus we have ^-= — ~~ whilst 

dn dr 


e_ n\ _ _d_ m _ j_ _ j_ 

dn \r/ dr \r/ r 2 e 2 


Further, dS = e 2 dQ, where dQ is the solid angle subtended by dS at 
O. Hence 

*KeK8— 

—> 47 t/o as s 0, 

where/ 0 denotes the value of/at O. Let e tend to zero; the result is 
Green’s identity 

--c£[4G) - x M dS ~r„Lv dV - 

If O is outside the region V contained by S t equation (56.3) 
immediately yields 

i[fi^) J r 8 £] dS+ S,V dV ^- 

Consider the case when the function / is harmonic, that is Af= 0. 
We deduce from (60.1) that 


7—V.A. 
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Accordingly, the values of a harmonic function / and its normal 

Plf 

derivative ^ at all points on the boundary of a closed region deter¬ 
mine the value of the function at all interior points. 

It can be shown that either a knowledge of/or ^ on S determines 

the other when / is harmonic. Consequently this result has little 
practical value. 

Further, consider the cases when S is a sphere with centre at 0. 


Since we have in virtue of (57.1) that 



That is, if a function is harmonic in a spherical region, its value at 
the centre is the arithmetic mean of its values on the surface of the 
sphere. It can be shown that a function with this property is harmonic. 

In Section 42 we saw that Gauss’s theorem can be extended to the 
infinite volume outside a closed surface provided that lim R 2 a.— 0. 

B-» oo 

In the same way Green’s theorems can be extended to the infinite 
volume provided that lim R 2 f grad g = lim R 2 g grad /= 0. Corre- 

B-* oo B-* oo 

sponding to g — ~ we have grad £ = — r and these conditions are 
lim ( —/r)/r = lim r grad /= 0. 

oo r-*oo 

Thus, under these conditions, we may apply Green’s identity to 
the case when V is the region outside the closed surface <S containing 
the point O. Now let shrink to the point O and we obtain 


4t r 



dV , 


where the volume integral is taken over the whole of space. 




CHAPTER 11 


Orthogonal Curvilinear 

Coordinates 

61 Curvilinear Coordinates 

Let us change from rectangular Cartesian coordinates x, y and z to 
curvilinear coordinates u, v and tv by means of the three equations 

x = f(u, v, to), y = g(u, v, tv), z = h(u, v, tv), 

where /, g and h are single-valued functions with continuous partial 
derivatives of the first order in some given region. The functions /, g 
and h must be independent and the condition for this is that the 
Jacobian determinant 

dx 
dtv 

By 
dtv 

8z 
dtv 

In practice this condition may not apply at some points and particular 
care may be necessary there. 

Under the given conditions it can be shown that u, v and tv can be 
obtained as single-valued functions of x, y and z with continuous 
partial derivatives of the first order. 

If one of the coordinates u, v and tv is held fixed and the other two 
allowed to vary, the point (x, y, z) traces out a surface called a 
coordinate surface. There are three families of such surfaces 
corresponding to constant values of u, v and tv. 

If two of the coordinates u, v and tv are held fixed and the third 
coordinate allowed to vary, the point (x, y, z) traces out a curve 

89 



dx dx 
8u dv 

d(x, y, z) ^ d)t dy 
d(u, v, tv) du dv 

dz dz 
du dv 
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called a coordinate curve. There are three families of such curves 
corresponding to varying w, v and to respectively. Each coordinate 
curve is the intersection of two coordinate surfaces. For example, the 
curve corresponding to varying u is the intersection of the surfaces 
corresponding to fixed v and to. 

Let P be given by the position vector r = xi +yj + #k where i, j and 
k form a fundamental system. The change of coordinates to w, v, to 

makes r a function of u , v and to. Then the vectors ^ are 

cu dv dto 

tangent to the coordinate curves. Denoting the unit vectors tangent 
to the coordinate curves by I, J and K, we have 

dr 

dv 


8r 

to - ***• 


— h 2 J, 




where 


h ! 


ho. — 


h* — 


dr 


dto\ 


62 Orthogonal Curvilinear Coordinates 

We shall now impose the restriction that the coordinate curves are 
mutually orthogonal. In this case we refer to an orthogonal 
curvilinear coordinate system and the vectors I, J and K satisfy the 
relations 


Further, we have 


J.K = K.I = I.J = 0. 


, dr , dr , dr , 
dr = — du + — dv + dto 
du dv dto 


= h x dul + h 2 dvj + h 3 dtoK 

and so the length ds of arc between the neighbouring points (w, v, to) 
and (u + duy v + dv, to+dto) is given by 

ds 2 = dr.dr = h 2 du 2 + h 2 dv 2 +h 2 dto 2 , 

since I, J, K form an orthogonal system of unit vectors. 

The arc lengths ds lt ds 2 and ds 3 along the coordinate curves are 
given by 

ds x = h x du, ds 2 = h 2 dv , ds 3 = h 3 dto. 

Example 1. Obtain h lt h 2 , h 3 and the vectors I, J, K for the 
cylindrical coordinate system given by x—u cos v, y — u sin v, z=to. 

The ranges of u, v and to are given by 0, 0 < z; < 277 and 
— co<to <oo. 
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11 ] 

We have 
Thus 


dr 

du 

dr 

dv 

dr 

duo 


r = u cos v i + u sin v j + zvk. 

\dr\ 


— = cos v i+sin v 1 , 
du J 


— = — u sin vi + u cos v i, 
dv 


= K 


\du\ 

dr\ 

dv 

dr 


dzo\ 


= 1 , 
= u, 

= 1. 


Hence 

and 


K = 1, h 2 = u, h 3 = 1 

I = cos v i + sin ©j, J = — sin v i + cos v j, K = k. 


Example 2. Obtain h ± , h 2 , h Q and the vectors I, J, K for the spherical 
polar coordinate system given by x=u sin v cos zv, y~u sin v sin zv, 
z=u cos v. 

Here the ranges of u, v and w are given by 0, and 

0 ^ zv < 2tt. We have 


r — u sin v cos zv i + u sin v sin zvj + u cos v k. 

= 1 , 


Thus 

dr . . . 

— = sin v cos zv i + sin v sin zvj + cos v k, 

dr 

— — u cos v cos zv i+u cos v sin zv i — u sin v k, 

dv J 

dr .... 

— = — u sin v sin zoi + u sin v cos zv i, 

dzo J 

Hence 


dzo 


= u. 


= u sin v. 


and 


= 1, h 2 = u, h 3 = u sin v 

I = sin v cos zv i + sin v sin zv j + cos v k, 
J = cos v cos zvi + cos v sin zv j — sin v k, 
K = — sin zv i + cos to j. 


Exercise 

1. Obtain h lf h 2t h 3 and the vectors I, J, K for the following coordinate 
systems: 
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(i) parabolic cylindrical coordinates given by 

x = £(u 2 —v 2 ), y = uv, z = tv, 

(ii) paraboloidal coordinates given by 

x = uv cos tv, y = uv sin tv, z = ^(« 2 — v 2 ), 

(iii) prolate ellipsoidal coordinates given by 

x = a\/{(u 2 — 1)(1 — ® 2 )} cos tv, y = ay/{(u 2 — 1)(1 — ® 2 )} sin tv, 
z — auv. 


63 Gradient 

In Section 38 it was proved that 

grad/= Jl + gj+^K. 

However, the directional derivative 3fjds 1 is given by 

= lim f(u+du, v, to)-f(u, v, to) = 
ds x 4u-»o h x du h x 3u 


c . .. . df 10/ ,0/10/ 

Similarly = — and = and 

USq /?2 OSg ^3 


SO 


era d/=i^I + i^T + i^K 

g d -' h x du 1+ h 2 dv i+ h 3 dtv i ^ 


(63.1) 


From this result we see that 


111 

grad u = j— I, grad v = j-J, grad to = j- K. 

«i h 2 «3 

Hence 

I = J x K = h 2 h 3 grad v x grad to,' 

J — K x I = h 3 h x grad to x grad u, * 

K = I x J = h x h 2 grad u x grad v. 


(63.2) 


Using the results of examples 1 and 2 of the previous section, we have 
that grad / in cylindrical and spherical polar coordinates is given by 


and 


0/ _ 1 0/ , 0/ Tr 
5-I+- j-J+t-K 
cu ucv oto 


0/ T 10/. 1 0/ Tr 

7r~ I H 7T" J H-:- ^— K 

cu u cv u sin v oto 


respectively. 
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Exercise 

1. Obtain grad/ in the following coordinate systems: (i) parabolic 
cylindrical coordinates, (ii) paraboloidal coordinates, (iii) prolate ellipsoidal 
coordinates. 


64 Divergence 

Consider the vector a= a t \+ 0 jJ + where the a It and a K 
are now components with respect to the basis I, J, K. From 

div(/a) = /div a + a.grad/ 
we have in virtue of (63.2) that 
div (< 2 /I) = div (ajh 2 h 3 grad v x grad to) 

= ajh 2 h 3 div (grad v x grad to) + grad v x grad to. grad (ajh 2 h 3 ). 
But div (a x b) = b . curl a — a. curl b and curl grad /= 0. Hence 
div (/Z/I) = grad v x grad to. grad (aji 2 k 3 ) 


1 


h 2 h 3 


J xK -[i 


1 8{a I h 2 h 3 ) 1 8{a I h 2 h 3 ) 


t ■ ~ t 

8u h 2 8v ^ 


1 8{q.Ji 2 h 3 ) 


8zo 


K 


1 ^ f 1 7 \ 


Adding the corresponding results for div («J) and div (a K K) we have 

diVa = hJhfh \^i ( aih2h ^ + ( a J h M KM 2 )] • (64.1) 

In particular, the corresponding results for cylindrical and spherical 
polar coordinates are 


and 


1 T 8 . x 8aj 8 , J 

[£ ( “ 2 ^ °°')+h (“ 8in ' s (■■«>] 


respectively. 

The Laplacian of / has been defined by Af= div grad/ and so by 
the combination of the formulae for gradient and divergence we have 

' f * » + i M . 1 M I\1 ,64,-, 

h x h 2 h 3 [du \ 8u) 8v \ h 2 8v) 8w\h 3 8w) J* ' ' 
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The corresponding results for cylindrical and spherical polar 
coordinates are 

u 8u \ 8uJ u 2 8v 2 8w 2 

and 

1LU2S \—*L 

u 2 8u \ duj u 2 sin v 8v \ 8vJ u 2 sin 2 v 8w 2 
respectively. 

Exercises 

1. If u t v, w are cylindrical coordinates, show that grad log u and grad v 
are solenoidal vectors in any region which does not contain the origin. 

2. Obtain div a and Af in the following coordinate systems: (i) parabolic 
cylindrical coordinates, (ii) paraboloidal coordinates, (iii) prolate ellipsoidal 
coordinates. 


65 Curl 

From curl (/a) =/ curl a + grad /x a we have 

curl («jl) = curl {a^h x grad u) 

= a^i x curl grad u + grad (ajhj) x grad u. 

But curl grad u = 0 and so we have 

curl (a,l) = {~ + 11 (*,*!»+£ ^ («A)K} x I 


1 8 

h 3 hx 8w 


(«A)J 


1 8 


h ± h 2 8v 


(«A)K. 


Adding the corresponding results for curl (tfj) and curl (tf^K) we 
have 

curl a = < aA) } J 
+ iS (aA) 4,^^} K - t 65 - 1 ) 

In particular, the corresponding results for cylindrical and spherical 
polar coordinates are 

1 r 8a K 8ajl T \8a t 8a K 1. ,1 \ 8 . 8aA 

U J+ U 


11 ] 

and 
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Exercise 

1. Obtain curl a in the following coordinate systems: (i) parabolic 
cylindrical coordinates, (ii) paraboloidal coordinates, (iii) prolate ellipsoidal 
coordinates. 


66 Curl of the Curl of a Vector 


With reference to a fundamental basis i, j and k we have 
. (da k daA. (da, da k \, (da* da,\ t 

cuda = {-^-^y + {ad-^P + {-d'^) k - 

and so the coefficient of i in curl curl a is 


d 

(da, 

da t \ 

d 

(da. 

dv 

l du 

dv J 

dzo 

\8w 


du) 

(d 2 a, d 2 a, d 2 aA d (da, da, 

\du 2 dv 2 duo 2 ) du \du + dv + 


daA 
dw J 


Let us define the Laplacian operating on a vector by 


Then we have 


Aa = Aad+Aafi+Aa^s.. 
curl curl a = grad div a— Aa. 


( 66 . 1 ) 


It is important to note that this relation has only been proved for 
a fundamental basis. Some authors define Aa in an orthogonal 
curvilinear coordinate system by this relation. The reader is referred 
to Section 76 for a full treatment of curl curl a. There we show that 
this relation is always valid for one particular basis, called the natural 
basis. 
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Miscellaneous Exercises 
1. Show that 


81 

8u 

2 

8u 

8K 

8u 

^Hint: The derivatives are linear combinations of the basis vectors. Use 
relations of the type 

a 2 r \ 


- 

1 8h x 
h 2 &v 


8_K „ 81 

8zv * 8v 

1 

~ h : 

8k 2 

L a« 

J, 

81 

8zv 

II 

8h 3 

8u 

K, 

1 

8h x 

3J = 

1_ dh 2 

1 

8h, K 


1 

8h 3 

K 

h 2 

8v J ’ 

8v 

h x 8u 

h 3 

dzv 

—7 

8w 

h 2 

8v 

AX, 

1 

8h x 

8K 

,±8h 2 

8K 


1 

8h 3 

1 

8h 3 


h 3 

7S * 

8v 

h 3 dzv 

dzv 

— 

K 

—-1- 
8u 

h 2 

8v 

J- 


I — =0 I — = 
8u * 8u 


*i“ d 


2. If I, J and K are the basis vectors for a spherical polar coordinate 
system, show that 


curl curl 1 = 0, curl curl J = cot ^ I ( 


curl curl K = -■ K, 


, _ 2 _ , j. T cot v _ cosec 2 v 

grad divl = —^ 2 1, grad div J =-^ 1 - 


J, 


grad div K = 0. 

3. In an orthogonal curvilinear coordinate system for which h x — l t h 2 is 
a function of u and h 3 is a function of u and v, show that 

i r 8 /h 2 h 3 ai\ , a (k 3 k x ai\, a (h x h 2 ai\i 


CHAPTER 12 


Contravariance and 
Covariance 


67 Contravariant Components 

Consider the basis formed at each point by the three independent 
vector fields e 1} e 2 and e 3 . These vectors are not in general unit vectors, 
nor are they mutually orthogonal. We may write 

a = a 1 ^ + a 2 e 2 + « 3 e 3 

and then refer to a 1 , a 2 , a 3 as the contravariant components of a 
with respect to the basis e a . For brevity, the convention is made that 
all Greek indices have the range 1, 2, 3. 

Further, we introduce the Einstein summation convention that 
a repeated index implies summation over the range of values 1, 2 and 
3. That is, we may write 

a = a a e a . 

Let the contravariant components of a with respect to a second 
basis e' be a' a . That is, 

a = 

Both sets of vector fields e a and e' a form bases and so we may express 
each vector of either basis linearly in terms of the vectors of the other 
basis. Accordingly, 

e ' a = pfet and e„ = (67.1) 

where p%, and q B a each represent nine quantities which are not in general 
constants. 

The repeated index j8, called a dummy index, in the above 
equation can be changed to y whilst the free index a can be changed 
to jS. As a result we may rewrite the equations in the forms 

e (t= Pl e v and e e = q$e' Y . 

97 
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Substitute in the previous equations to obtain 

< = and e a = p}q e a e vy 

where the two sets of repeated indices imply double summations. 
Introduce the Kronecker delta by 



_ fl if y = a. 


“ \0 if 

and we deducef that 

ptal = ptf. = K- 

Further, we have 

a = a a e a = a' a e,' ay 


and so 


a a qie’(i = a ,p e' p and a*e 0 = a ,cc p^e e . 

The vectors e # and e# form bases and so no linear combination of 
them can vanish. Accordingly, we have 

a' 0 — qlfP and oP = (67.2) 

68 Covariant Components 

Let the vector fields e 1 , e 2 and e 3 be reciprocal to e lt e 2 and e 3 . 
Then the conditions of (13.1) can be written 

e a .e p = 8£. 

If e'“ form the basis reciprocal to e^ we have 

e a-e //J = 8 ' 

We may write 

a = Gje 1 + a^e 2 + a 3 e 3 = a a e a . 

The quantities a u a 2 , a Q are called the covariant components of a 
with respect to the basis e a . Note carefully that a ly a 2 , a 3 are con- 
travariant components with respect to the reciprocal basis e“. 

From a = a a e a = a! a <e! a , we have in virtue of equation (67.1) that 

a a e a .e p = a' a e' a .e p = <e' a .$e; 

and so 

= a' a q}S%, 

from which by summation we have 

a p = ffX- (68.1) 

t The reader familiar with matrices will recognize that [q%\ is the matrix 
inverse to M- 
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Similarly, we obtain 

a'p = ppCiy. ( 68 . 2 ) 

Further, we have in virtue of equation (67.2) that 


a = a Y e r = a' e e' p = p%a Y e' e 

and so 

aW-pie") = 0 . 


The components a Y are arbitrary. We may in turn select one com¬ 
ponent to be unity and the others to be zero and we obtain 


Similarly, we obtain 


e y = ple ,e . 
e' y = 


(68.3) 

(68.4) 


Exercises 

1. Show that a B =pla Y . 

2. Show that e' Y =q}e p . 


69 Fundamental Tensors 

We define g a p=gp a * called the covariant fundamental tensor, by 



Saff 

— ©a • ©0 gfiw 

From 




a 

= a a e a = a a e a 

we have 




o“e a .e p = 

a a e a .e fi = a a h\ 

and so 




«*= gafi<P- (69.1) 

Further, define g aB = g ea , called the contravariant fundamental 
tensor, by 

g“t = e a .e e = g Ba . 

Then we have 


a“e c .e^ = a a e a .e B , 

from which we obtain 


a? = g° e a a . (69.2) 

Again, using equation (69.1) we have 

a a e a = a B e ff = g afi a a e e 

and so 


8 —V.A. 


a a (e a -g ae e e ) = 0. 
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Since the components a a are arbitrary, we have as in the preceding 
section that 


= g a0 e B . 

In a similar way, we obtain 

(69.3) 

e" = g“ B e 0 . 

Hence we have 

(69.4) 

e“ = g aB g Bv e y 

and sof 


g a %y = 

(69.5) 

Exercises 



1. Show that e a =g aB e e . 

2. If a and b are vectors, show with the usual notation that 

a.b = d*b a = aj> a = g aB a u b e = g aB a a b B . 


70 Natural Basis 

Consider the coordinate systems u a and u a and let r denote the 
vector position relative to the origin O of the u a system and r + c 

denote the vector position relative to the origin O' of the u' a system, 

— > 

where c is the constant vector O'O. We choose as the natural bases 
the vectors e a and e' a given by 


= 


dr_ 

du a 


and 


, _ d(r + c) dr 
a “ 8u' a ~ du’ a ' 


The chain rule of partial differentiation states that 

dr _ dr du B dr _ dr du' B 

du >a du B du' a an du a du' B du a 

Hence we have 


e' = 


du B 

dH“ 


and 


— 


du' B , 
du a 


Accordingly the p B and q B defined in Section 67 are given by 


P 


B _ 
a 


du B 

du' a 


A B ^ U ' B 


f That is, the matrix [g a/? ] is the inverse of the matrix [g a/3 ]. 
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and so, using (67.2), (68.1) and (68.2), the contravariant and covariant 
components of the vector a transform by 


a p = 


a p = 


— a'* 
0 m '“ ’ 

'B du'* 
a ~ 8u a a ’ 

8u' a , 

, 8u tt 


ap ~ 8iP a “* 


(70.1) 


The transformation laws with respect to the natural bases have a 
simple structure and provide the starting point in the development of 
the tensor calculus. 


71 Physical Components of a Vector 

We have defined the contravariant and covariant components of a 
vector with respect to any given basis. But the component as 
defined in Chapter 1, for emphasis now called the physical compo¬ 
nent, of a vector a along the direction given by the unit vector 1 is 
a.l. 

We have e a .e /5 =^ a/5 and e tt .e p —g aP and so |e a | 2 =£ aa and 
\e tt \ 2 =g aa , where there is no summation over the repeated index a. 
That is, the unit vectors in the directions e ff and e a are respectively 
e JVg«« an d e a IVg aa > ( no summation). Hence the physical compo¬ 
nents of a along the directions e a and e a are 

a • ^al\^Saa = • ®a/'Y/<§ , aa = = ^a/Vj§aa 

and 

a.e«/V£ aa = «%.e a /V£ aa = a p h^g aa = a a Wg aa . 

When the basis vectors consist of three mutually orthogonal unit 
vectors we have g aP — 8% and consequently e a =e®. That is, the basis 
is self-reciprocal. We see in this case that the contravariant, covariant 
and physical components are identical. 

Consider the vector grad/ whose components referred to a 
rectangular Cartesian coordinate system x a are df[dx a . On trans¬ 
formation to the u a coordinate system, the covariant components with 
respect to the natural basis are given in virtue of (70.1) by 

8f dx? 8f 
dx a 8u B du e 


If the u a form an orthogonal curvilinear coordinate system, we 
have g aa —h 2 (no summation) and so the physical components of the 
1 

gradient are (no summation) as already obtained in Section 63. 
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72 Derivatives of Natural Basis Vectors 

The natural basis vectors have been defined by e a =dr /8u a and so it 
follows that 

_ 8e 0 

8uP 8u a 

Further, differentiation of the relation ^ a .e fi =g afi with respect to u y 
yields 

tea _ , . _ dg a3 

8u Y ^ “ 8u Y 8u v * 

Note carefully that 8eJ8u Y represents nine vectors, each of which can 
be expressed in terms of the basis vectors themselves. Accordingly, 
we may write 

tea — X p e 
8vP ~ p ’ 


where represents the 27 coefficients involved. Substitution in the 
previous equation yields 

"b e a • = ^y * 

That is, in virtue of the definitions in Section 69, 

gpfi^ay ^~gap^$y = ~jfay' 

Permuting the indices a, jS and y cyclically we obtain 

gpytfia +g0p^ya = 

and 

gpc^fo+gypKe = f^T 


We introduce the Christoffel symbols [ajS, y] and 



defined by 


and 


(dgfir 

1 dgay 

Qgafi\ 

\8u a 

+ 8u* 

8u Y / 




In virtue of 


8e a _8e fi 
8u* 8u a 


we have X P B = X p Pa and so a simple calculation 


12 ] 

gives 
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Spv^afi —■ [ a A y\' 

Hence we have 


K, = KK, = g”g,X> = y] = {^}. 

and so the derivatives of the basis vectors are given by 


8u B 



(72.1) 


To obtain the corresponding derivatives of the e c , we differentiate 
the relation e a .e B = S B with respect to u v to yield 


e« 


8e B 8e a 
8u? + 8tf 


• e* = 0. 


8e B 

Let ■—-=u^„e x> and we have by substitution that 
8u y 


e„.e» = 0. 

That is, 

«» +s »ty} = ° 

from which 



Accordingly the required relations are 



(72.2) 


73 Derivatives of Vectors 

Consider the vector a= a a e a . On differentiation we have 

da 8a a „ 8e a 
8u B 8u B a+ 8u B 

8a° f p 'l 

"5? e » + *ta; e » 

= a’jfin 
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where we have written 



The expression a p 0 is called the covariant derivative of a p with 
respect to u 0 . 

Similarly from a = a a e a we have 


da da a de a 

du* du* + a du* 


where we have written 



a 



e p 


= a 


P.0* 


a p.B ~ 


da p f a\ 
du e \ppf 


a 


a* 


The expression a pS is called the covariant derivative of a p with 
respect to u e . 


Exercises 


1. If the basis vectors are constant, show that covariant differentiation is 
partial differentiation. 

C^2l C^2l 

2. By considering ■ = ■ . show that 

& 8u a du e 8u e du? 


d_ 

du v 




= 0. 


74 Gradient 

In Section 71, we saw that the covariant components of grad / 
referred to the natural basis e a are df/du a . Accordingly, we have 


grad/= ^ e “- 


(74.1) 


Hence 

grad u a = e“ and curl e“ = curl grad u a = 0. (74.2) 

Let us write A = [e 1 e 2 e 3 ] and we have 

I?-[£ eV ]+[ ei £ e3 ]+[ e ' e2 £] 

= -Q M-ft} [ eW Hi} [ eW ] 



12] 

Consequently, 

we obtain 
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grad A = 

II 

a 

R 

II 

1 

i} Ae " 

(74.3) 


75 Divergence 

From equation(46.1) we have div (e“ x e^) = . curl e“—e“.curl e p . 

Further, equation (74.2) states that curl e“=0 and equations (13.3) 
take the form Ae x = e 2 x e 3 , etc. Accordingly we derive div (Ae a ) = 0 
from which in virtue of (43.1) we have 

A div e a +grad A.e a = 0. 

On application of (74.3) we deduce that 

Using equation (69.4) the divergence of the reciprocal vectors may be 
calculated from 


div e“ = div (e“.e^ e p ) 

= e a -e p div e fi + grad (e“. e^). e p 

^0® r cc 

Use of e y . = 8|, = — -l > e a and appropriate changes of dummy 

indices yields 

dive«= -**{£}. 

Accordingly, we have 

div a = div (a a e a ) = div (a a e a ) 

= a a div e a + grad a a .e a = a a div e“+grad a a .e“ 
da a . 


[pep) 

$0 4 -1 @ \ n a 

a + \«P) 


e p .e« 


= ~ a ag P7 


da a 

du* 


g' 


U da a 
+ fa* 

(&-&W 


e^.e“ 


— n tt 

u ,a 


l a.fi 
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We deduce that 


[Ch. 


grad div a = — (a a>0 e a .e e )e Y 



= a a l3y e a .e /J e y , 

where we define a a By by the expression in the square brackets in the 
line above. 

The Laplacian da of the vector a is defined by 
da = « a>j8y e a e*.e y = g eY a a ^ y e« 

The reader is asked to verify that this definition produces the same 
result as that of Section 66 in the case of a fundamental basis. 

Hence we have 


grad div a—da = a a>ey e a .e e e Y —a aJ}y e a e^.e 7 
= « a jgy x (e y x e a ). 


Exercises 
1. Show that 


Aa=g 


( 


g 2 a 

8u p du v 



2. Using the result of Exercise 1 of Section 73, show that a aeY =a aYB . 


76 Curl 

Since curl e a =0, we have in virtue of (45.1) that 
curl a = curl (a a e a ) 

= a a curl e a +grad a a x e a 


_ 


du Y 


e 7 x e“. 


Since fc} is symmetric in a and y and e y x e a is skew-symmetric in 
a and y, it follows that < ° >e y xe“=0 and so 

\ccyj 
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Further, 


curl curl a = curl (« a y e y x e a ) 

= curl (a a Y e v x [e a .e ff ]e ff )f 

- {^K, e 'x e °- e «)} e "x e ” 

- {ihs [ e ’ e ° e J J £ }t e ' e “ eJ }[ e ' e ' e J 

+ °-u} [e,e “ eJ } e ' xe " 

The last summation vanishes since -I e >e fl xe (7 =0. Hence 

\?P) 

curl curl a = [e’e“e„] {^- {^} ««.» - {£} e» x e» 

= 0 a>w >[e y e a e ff ]e* X e ° 

= a a.yp eP x {e y x [e a . e ff ]e ff } 

= a a . Y pe p x (e y x e c ).f 

In virtue of Exercise 2 of the preceding section a awVP is symmetrical 
in y and p. Hence we have 

curl curl a = a a py e p x (e y x e a ) 

= a a l3y e e x (e y x e a ). 

That is, we have established that 


curl curl a = grad div a—Ja. 

It is well to emphasize that this relation holds for the tensor com¬ 
ponents referred to a natural basis, but not in general for the physical 
components of the vector a. 


t(e“.e a )e* = S«e* = e“. 




Solutions 


Chapter 1 

§ 1 1. (c) and (/) are vectors, others are scalars. 

§ 6 1. (i) b—a, (ii) —a, (iii) — b, (iv) a—b, (v) —2b, (vi) 2a—2b. 

Components of FA are — 5\/3/2 and 5/2. 

Components of FE are 5\/3/2 and 5/2. 

2. AB= a+P+Y> CA= —a—P, projected vector is — (a+P)/2. 

§ 7 1. (i) 6i—3j + 2k, (ii) 3^/5, (iii) -(6i+6j-5k)/V(97). 

§ 8 2. (ii) sin" 1 (31/50), (iii) ± (i — 11 j — 7k)/V( 

§ 9 2. (i) — 3(i+j+k), — (2i+j+2k), -(i+2j+4k),(ii) ±(2i+j + 2k)/3, 

(iii) V(21)/2. 

§10 1. 1. 

§ 11 3. (i) -22, (ii)27i+4j-Ilk, (iii) -2(5i+j-9k), (iv)ll(2i+4j-6k) 
(v) 540, (vi) -1 l(54i+8j -1 lk). 

§13 1. Mj+k),i(k+i),i(i+j). 

Miscellaneous Exercises 
1. (bxc)/[abc]. 

Chapter 2 

Miscellaneous Exercises 

1. 7. 

2. 9/V(17). 

3. 58/(5V2). 

9. (a) cos -1 (b) cos -1 1 ly/3. 

Chapter 3 

§ 19 2. p=£tt 8 a+itt 2 b. 

§ 23 2. (i) (1 + u 2 ) sin u—u cos u, 

(ii) (m 8 sin u—Zu 2 cos «)i+(tt 8 cos u+3u 2 sin u) j 
+ (cos u+u sin u+u 2 cos «)k. 

§ 24 1. (i) ve uv i + 2j+ v cos u k, (ii) v 2 e uv i — v sin u k, 

(iii) (1 + uv)e uv i + cos u k, (iv) u 2 e uv i, 

(v) (2 sin u+v cos u) i+ v(u cos u —sin u)e uv j — (2 u + v)e uv k. 

2. (i) 0, (ii) 
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Chapter 4 

§28 1 . (—/cos«?i—/sin «?j+k)/\/(l+/ 2 ). 

§ 29 1. a 2 (J« 2 +sin 2 « dv 2 ), (1 +f 2 ) du 2 +u 2 <fo 2 , rf« 2 +-|« 2 + j- <fo 2 . 

3. — a( J « 2 + sin 2 © <fo 2 ), (/ du 2 + uf dv 2 )[( 1 +/ 2 ) 1/2 , 


2 )/(“ 2 +(l) 2 ) 


2 \ 1/2 


Chapter 5 

§32 1. (i) 0, (ii) 0. 

2. (i) 2, (ii) 3, (iii) 2flfe 
§33 1. p(c 2 *-l)k. 

§35 1. 6(i+j+k),0,O. 

2 . 1 / 8 . 

§36 1. 24i+96j+394k/5. 


Chapter 6 

§ 37 1. (i) 7^6/2, (ii) -9/V2, (iii) lIV 2 - 2. 2A.o. 
§38 1. (i) 2, (ii) V3, (iii) 

§ 39 3. No, (ii) Yes, £(X.p)r 2 , (iii) Yes, (A.r)(p.r), (iv) No. 


Chapter 7 

§ 41 3. (i) 2(x+y+z), (ii) 0, (iii) zx, (iv) 0, 

(v) yz(2x +y + z) i+ zx(x +2y+z) j -f- xy{x +y + 2z)k, 

(vi) 0, (vii) 2zx(z—x)(x 2 +6y 2 +z 2 +zx). 

4. (i) 0, (ii) 6, (iii)^ + ?|,(iv)l 

§42 1. 0. 

3. 47r(a+6+c)/3. 

§43 1. cjr 3 for any constant c. 

Chapter 8 

§ 44 1. (i) 1+j-k, (ii) -i+j-t, (iu) 0, (iv) 0, 

(v) (3 y 2 z+3yz 2 +2xy 2 +2xz 2 )i+z(3x 2 + 3 zx—2yz)j 
-\-y($x? 4- 3xy—2yz)k, (vi) 0, (vii) 0. 

§45 2. (i) e t+v * !t {x{y—z)i+y(z—x)l+z{x—y)k)> 

(ii) x(ze* —ye v )i +y(xe x — ze*)j + z(ye v — xe x )k. 


Chapter 9 

§ 50 1. See equation (44.3). 

§ 55 2. (i) (y 2 + z 2 )i + (z 2 + ac 2 )j +(# 2 +y 2 )k, (ii) r 3 r x (r x c)/6, 

(iii) %{x(z 2 —y 2 )i +y(oP—z 2 ) j+ z(y 2 — x 2 )k), 

(iv) (xe x -f- 1 — e x )(zi—xk)/x 2 . 
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SOLUTIONS 


Chapter 11 

§ 62 1. (i) h x =h 2 = V(« 2 +^3=1, 

I=(ai+®j)/\/ (u 2 + ® 2 ), J=(—®i+ «j)/V(« 2 + v2 )> K=k, 

(ii) h 1 =h 2 = '\/ (a 2 + v 2 ), h 3 =uv; 

I=(® cos w i+v sin a; j + uk)/^(u 2 + v 2 ), 

J=(a cos zv i+a sin a?j —®k)/\/(a 2 +® 2 ), 

K= —sin zv i+cos zv j. 

(iii) h i = «V((« 2 -® 2 )/(« 2 -1)}, K = a^/{{u 2 -v 2 )l(\-v 2 )} y 
h 3 =«V((« 2 -1)(1 - v 2 ) ; 

i= “V(^$) {cos ” i+sin 
j=K,i+sin 
K= —sin a? i+cos a? j. 

*« *• (I i+ Ij)m“ 2 +^)+Ik, 

(iii) |v(“ 2 -i) I I+V(l-^ 2 ) |j}/W(* 2 -*- 2 » 

+ gK[^ 2 -lXl- 2 )}]. 

§ 64 2 - <9 {|;(v'(“ 2 +®> 1 ) +|, (V(“ 2 +e a K)+(-< 2 +® 2 ) 

i fa 2 / ay\ ay 

a 2 + v 2 \du 2 a ® 2 J dzv 2 * 

(ii) uv(J+v=) {m [“®v'(“ 2 +® 2 )«/]+|; [«V(“ 2 +*>>,]} 

_j_1_ & a K 

a® 8zv 

i r_a / an a/ 8f\\ 1 ay 

uz)(u 2 +a 2 ) \a« \ D a«/ a® \ u a®/ j a 2 ® 2 aa? 2 ’ 

< m ) {I (V«« 2 -iX- 2 - 2 )K) 

+(V{( 1 -® 2 )(" 2 -'’ 2 )} «/)| 

1 

«\/{( w2 — i)(i—® 2 )} aa? ’ 

, 1 g 2 / 

« 2 (a 2 —1)(1 — ® 2 ) dzv 2 
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§65 1. (i) -!_ . fk? 

\Z(u 2 +v 2 )\ dv du j dzt 


^ uv-^(u 2 +v 2 ) 8v a 
uv \dw * 8w j 


(V (« 2 


J 


+ 

uv 

1 

^h 2 +w 2 


w vw j 

\ {J; (V(« 2 +® 2 K)-^ (V(« 2 +® 2 K)} k, 
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VECTOR ANALYSIS 


This revised second edition provides a systematic accoun o ^ 
algebra and calculus from the fundamental definitions. Improvements 
include a rigorous proof of Taylor’s theorem for a vector function, a 
vectorial definition of directional derivative, more material on Green s 
identity and several extra examples. The book is intended for first year 
students reading an honours course in mathematics and for students ot 
mathematics, the physical sciences and engineering at universities and 
technical colleges. 

The theory is developed with emphasis on the vector as an entity in 
itself rather than on its definition as a number triple. It includes a chapter 
on line integrals, surface integrals and volume integrals, as many students 
require these concepts in vector theory before they have taken a lorirnl 
course on integration. 

The theory is presented independently o f any co-ordinate system. 
In particular, the gradient is defined in terms of the directional derivative 
whilst the divergence and curl are defined by means o limits ot integrals. 
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